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Summary 


It is suggested that the heat lost by the corona is made good by energy 
given up in the viscous damping of hydromagnetic waves originating in 
the chromosphere. An estimate of the energy which the waves, generated 
by turbulent motion, would carry gives a value of 10° ergs cm~* sec™!. A 
frequency for the waves of 10~? sec~!, which seems to be the lowest possible, 
gives a fall in amplitude by 1/e in 10!° cm at a temperature of 10°°K. At 
higher frequencies, the distance for the same fall would be less. 

The damping is found to be less the smaller the inclination of the 
wave-normal to the magnetic field. For a dipole type solar field, less 
energy would therefore be given to the corona over the poles than over 
the equator. 





Introduction.—That turbulent motion in the photosphere would give rise 
to hydromagnetic waves, which would travel upwards into the chromosphere 
and corona and there give up their energy, was suggested first by Alfvén (z). 
This theory of coronal heating, in its original form, supposed that enough energy 
would be taken from the waves by Joule damping. It was therefore unable to 
withstand the criticism of Cowling (2), who showed that heating by electric 
currents in the waves must be far too small to maintain the temperature of the 
corona. It has been revived by the discovery by Piddington (3) that, under 
certain conditions, collisions between heavy ions moving with the wave and 
neutral particles can dissipate the energy of the wave much more quickly than 
Joule damping. In this case, the relative velocity of the colliding particles is 
not an electric current, bounded by limits to the strength of the magnetic field. 

Piddington showed that hydro magnetic waves would reach the chromosphere 
and there be damped, in the neighbourhood of transition between neutral and 
ionized hydrogen, strongly enough for the heat lost from all higher levels to be 
made good. Although his process in this way explains the transfer of heat against 
the temperature gradient from the photosphere to as far as the chromosphere, it 
cannot be carried higher because of the lack of neutral particles in the corona. 
As a step towards explaining the high temperature of the corona, he has pointed 
out that protons in the undamped wave in the chromosphere could have speeds 
~10’cm/sec and, in collision with neutral particles, these might be thrown 
upwards with an equivalent thermal energy to form a gas at 10° °K. 

An objection to this idea of upholding the corona comes from the short mean 
free-path of protons at a temperature of 10° °K and a density of 10®cm~*, which 
is about 10’cm, compared with the distance of at least 10°cm over which the 
temperature rises from the level where hydrogen becomes ionized to where 
10° deg. is reached. The temperature would reach its maximum at the point 
where the upward-moving protons undergo their first collision. 


*Received in original form 1960 May 6. 
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Some other way of transferring heat from the chromosphere to the corona 
seems needed. The dissipation of the energy of a wave by proton-electron 
collisions, as Joule damping, is known to be much too small. There remain 
proton-proton collisions, which lead to a loss of energy by the wave when the two 
particles are brought together by their thermal motion from different points in a 
non-uniform velocity field in the wave. This effect may be looked on as viscous 
damping. Although it is known to be negligible in the chromosphere, the 
coefficient of viscosity for collisions of this type increases as T°”. It is put forward 
in the present paper as a means by which hydromagnetic waves can heat the 
corona. 

Hydromagnetic waves with viscous damping.—In a gas made up of equal numbers 
of protons and electrons, each with number-density nm and partial pressure p, 
in which the particles gyrate much faster than they collide, the equation of motion 
of the protons may be written 


d e I 

ma, V= —VaH+cE——Vp+F. (1) 
Here V is the velocity of the centre of gyration, E and H the electric and magnetic 
fields and F the forces of viscosity and friction, taken to be small compared with the 
other terms. ‘There is a similar equation for the electrons. Neglecting F and 
supposing p to be given as a function of m, the operator V a used on this equation 


and the relation Va E= — Hf give 
9 
5 (= Vav+H) =VaVa("VAV+ Hl), 
ot\ e e 


Here the vorticity mc/e VaV would appear later as a Hall field, which is known to 
make the gas doubly-refracting for hydromagnetic waves. It is negligible 
compared with H for waves of frequency small compared with that of gyration 
of the ions. In this case 


<H- VaAVaH. (2) 
If A is a vector such that 


ViA=H 

this equation is unchanged by a gradient transformation of A, which may therefore 
be chosen so that equation (2) becomes 

é 

ot 
In general, V .A#o, since the gradient vectors in equation (1) which vanish under 
the operation V «are included in (0/dt)A and their divergence does not always 
vanish. 


When the momentum-density of the electrons is neglected, the equation 
of conservation of momentum for the gas is 


A=ViH=Vx(VaA). (3) 


d 
in which = is the Maxwell tensor, the electric field being neglected in the second 
form on the right-hand side, and II the pressure, made up of the isotropic pressures 
and the viscous tensors. Friction, or its equivalent resistivity, does not appear 


mms, V=V. (+I) =~ 2 HaVaH—2Vp+V.P (4) 
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explicitly in this equation. It falls out altogether when the momentum of the 
electrons is neglected. 
Suppose now that the magnetic field has the form 


and H,, to the first order, constant. A vector B can then be found so that 


Ia 
ey ae A= —-HyB. (5) 
Put into equation (4), these give a first-order linear fon al 
2 
mo B= a — Hy r{Va(VaH,aB)}+ 2ykTVV.B+ [ #v. he x8 |: 


Viscosity is here represented by a typical term only, put in brackets. Without 
terms of this type, the equation has solutions which are undamped waves moving, as 
Piddington (4) has shown, in a triply-refracting medium. This may be seen at 
once from the equation here as, if B is a plane wave of frequency », it is a set of three 
scalar equations homogeneous in the initial values of the components of B, in 
each of which there is a coefficient v?, giving in general three values of v? for a wave 
of given length and direction. 

As will be seen later, the waves chiefly of interest here are those for which 

H,2/421m, > ykT 
so that their potential energy lies in the magnetic field. Damping of waves of this 
type has been considered by van de Hulst (5) but without taking into account the 
effect of the magnetic field on the viscous stresses. This effect is not small and 
cannot be neglected, except for his purpose of showing the order of magnitude of 
the viscous term. 

Taking a rectangular system of co-ordinates, with H, along the z-axis, B may 
be taken to have a single component B(x, z, t) along the x-axis. A complete 
expression for the stress tensor in a magnetic field has been given by Chapman and 
Cowling (6), and the components which are at present relevant are 


a 7 tems et eee Vath, 
(+2 wt) 3 dx 6 dx 9 
9 
ee 2p 0 


=— ———-~ V 
zz xz 
(«+ 4 utr) dz 


Here w is the frequency of gyration of the particles and + the time between their 
collisions. The coefficient of viscosity for inverse-square collisions is 


w= 3 (=) (Gr) | Ai) aTS®, a=1x 107% 
7 





this value of «, which varies very slowly with the density, holding for both the 
chromosphere and the corona. As yp is proportional to m!?, the viscosity of the 
electrons is negligible. 

When the frequency of collision is small compared with that of gyration, P,, 
and P,,, become 


a 
Pi_= — | V, Pyg=0. 
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When these are compared with their forms for no magnetic field (w =o), it may be 
seen that the magnetic field has suppressed the force due to a shear in V, along the 
z-axis, which is remarkable because the free-path in this direction is unaffected. 
The effect is due to a tendency for ions in thermal motion along the field to be 
brought by electromagnetic forces to the local velocity of drift normal to the field. 
Putting V,, in P,,, in terms of B, equation (6) has the one component 

a HH? (2 o p & a 

Fy ar 47MXp (338+ 333) " waa (7) 
with solutions in the form of plane waves 





s 2 
B= By.exp.2mi (vt i. - iy + ies) ‘ 7 + % = = = a =5) 
a5 8 pv?x8 ae 2 mvs 
co 3 uH,? 3 mnu>’ 
If the wave-normal makes an angle @ with the magnetic field, the damping factor 
may be written 





Zz 
n€5- 


exp( — 27«z tan @sin @) 


so that a wave is damped most strongly when its normal lies along the x-axis and 
undamped when its normal lies along the z-axis. Equation (6) may also be put 
into cylindrical polar coordinates and solutions found as spherical waves radiating 
from the origin 

32)i2 


.[(p?+2 . 
B=B,| expani{ P22)" + vt+ ——_, Ero /@ +22)U2, 


The damping factor, when @ is the angle Bex wave-normal and the polar 
axis, has the same form as above, so that the amplitude is unchanged along the 
magnetic field and falls along any line at an angle to it. 

Transfer of energy from the chromosphere to the corona.—Measurements of the 
profiles of lines on spectrograms of the emission of the chromosphere have been 
made by Redman and Suemoto (7) and interpreted as showing line-of-sight velo- 
cities of the emitting atoms rising to 16 km/sec at 2600 km above the solar surface. 
By means of interferometric methods, Clube (8) has concluded that a turbulent 
velocity of 16km/sec is to be found at levels below 3000km. The level where 
hydrogen becomes ionized is between 6000 and 7oookm and, according to the 
theory of Piddington mentioned earlier, it is here that energy is taken up most 
rapidly from hydromagnetic waves. A mean speed of turbulent motion of 
20 km/sec therefore seems quite likely at :his level. 

Taking the density there to be 3 x 10°cm~* and the speed of hydromagnetic 
waves I x 10’cmsec~!, the energy carried upwards by waves arising from this 
turbulence would be 


1°7 x 1074 
3 x 10° 4+———_ (2 x 10°)? 1 x 10’= 1 x 10° ergscm~*sec™. 
2 
This is of the same order as estimates of the rate of loss of energy by the corona. 


1/2 
At this level T'~ 104 °K, so that (—) ~ 10®cmsec~}, 


Damping of the waves depends on their frequency. While there is no direct 
information about the period of turbulent motion in the chromosphere, an upper 
limit to the wavelengths formed there is set by the variation in density. If the 
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density m) in equation (7) is taken to vary exponentially along the z-axis 
with a scale-height , solutions of this equation, neglecting viscosity, which 
asymptotically become plane standing waves, are Bessel’s functions of the first and 
second kinds of order zero 


B= By, e**™ J, (s) and Bye Y,(s), s=4n€/A. 


When s <1, neither function is oscillatory and Y,(s)—> 0 asso. A wave-train 
can only be propagated if s>1, that is A<47¢. As in the chromosphere 


€=1'1 x 108cm, 
waves could not be formed to travel upwards with length greater than ~ 10°cm, 


corresponding to a frequency of 10~*sec~!. This gives the lowest value for the 
damping coefficient in the corona to be 


Kmin = 10-* # = 10- cm. 

As the speed of hydromagnetic waves is proportional to Hy, the upward flow 
of energy from the chromosphere will be greatest in regions where the magnetic 
field is unusually strong. For waves moving radially outwards, damping is more 
effective the greater the inclination of the magnetic field-lines to the vertical. If 
the steady field of the Sun is like that of a dipole, more energy should be absorbed 
by the corona over the equator than over the poles. A torsional oscillation during 
the sunspot cycle, by which the field-lines in the solar atmosphere are drawn out 
longitudinally to their greatest length at sunspot maxima, would at these times 
give the least difference between the polar and equatorial regions of the corona. 

Acknowledgment.—The suggestions of a referee of the Society have greatly 
improved the argument in this paper. 


50 Hans Place, 
London, S.W.1: 


1960 September. 
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FIRST RESULTS OF RADIO STAR OBSERVATIONS USING THE 
METHOD OF APERTURE SYNTHESIS 


P. F. Scott, M. Ryle and A. Hewish 


(Received 1960 September 20) 


Summary 


Four surveys, centred on declinations of 52°, 50°, 42° and 05°, have been 
made at a wave-length of 1-7 m with a large interferometric radio telescope 
employing the new technique of aperture synthesis. Details of the observa- 
tional method, calibration, data-reduction and analysis of the computed 
results are given. An assessment of the performance of the system shows 
that accurate measurements of the positions and flux densities of radio stars 
may be extended to limits well beyond those of any previous survey. Lists 
of sources are given for declinations between 52° and 50° and between 40° 
and 44°. Observations at o5° were made to study the distribution of radio 
stars outside the galactic plane; no evidence is found for a population of 
radio stars in the galactic halo. 





1. Introduction—The method of synthesizing a large radio telescope by 
combining observations made with movable elements in a number of different 
relative positions has already been described (Ryle and Hewish 1960). It was 
there shown how both pencil beam instruments and interferometers could be 
developed and how, in each case, one convenient system employed a long fixed 
aerial running East/West with a smaller aerial arranged to move on rails running 
North/South. With this arrangement it is only necessary to provide aerial 
rotation about the long axis to scan in declination while making use of the Earth’s 
rotation to scan in right ascension. 

A large interferometer operating at a frequency of 178 Mc/s which uses this 
method has been built for the study of radio stars. It was designed to obtain 
better positional accuracy for weaker sources so that a larger number might be 
available for optical examination and to obtain more reliable information on the 
apparent distribution of sources in space. ‘The same radio telescope has already 
been used in an investigation of the surface brightness of radio stars (Leslie 
1961) and for the measurement of the positions of 64 intense ones (Elsmore, 
Ryle and Leslie 1959). 

The present paper is concerned with the use of the instrument with its full 
resolution. The envelope of the interference pattern then has half-widths of 
25’ and 35’ in right ascension and declination respectively. The operation of 
the instrument and the methods of computation and data reduction are described 
in Section 2. The performance of the system is discussed in Section 3, and in 
Section 4 preliminary results of the observation of four areas of sky are summarized. 
The first two areas observed were centred on declinations of 50° and 52°; the 
later areas, observed with an improved recording system, were centred on 05° 
and 42°. 
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2. The radio telescope. 

(a) The aerial system.—The aerial system consists of two elements, one of 
which is fixed and the other mounted on railway lines. Structural details have 
been discussed elsewhere (Ryle 1960). In the present observations the aerials 
are used together as a transit interferometer. Both elements are cylindrical 
reflectors having a parabolic section 65 ft in width and each is aligned on an 
East/West axis. The fixed element is 1450 ft in length and the movable one, 
190 ft in length, traverses a distance of 1000 ft along rails in a North/South 
direction as shown in Fig. 1. 
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Fic. 1.—Schematic diagram of the aerials. 


Both elements may be set to declinations between —05° and go°. The 
reception pattern of each has a width of 4°-6 between half-power points in 
declination and this determines the area of the strip of sky under observation at a 
particular setting. In right ascension the reception pattern is related to the 
product of the reception patterns of the two aerials and to the frequency response 
of the output filter of the receiver as shown later. Its shape can be regulated 
either by the distribution of excitation along the fixed aerial or by means of the 
output filter response to give the narrowest reception pattern consistent with the 
lowest realizable side-lobe level. Side-lobes in right ascension at angles greater 
than 5° from the main beam are in any case removed by the lack of coherence 
between the signals arriving at the widely spaced aerials when using a receiver of 
bandwidth 4 Mc/s. 

Initially, the desired reception pattern in right ascension was achieved by 
adjusting the excitation of the fixed aerial. Since the spatial frequency spectrum 
of the reception pattern is derived from the convolution of the aperture distribution 
of the two aerials it was only necessary to adjust the excitation in steps of length 
equal to that of the movable aerial. The adjustment was performed by the use 
of unsymmetrical transformers at the junctions of the branched feeder network. 
The aperture distribution adopted, the effective spectrum of the combined aerial 
system and the resulting reception pattern are shown in Fig. 2. 

Later it was found more convenient to adopt a uniform excitation along the 
fixed aerial and to achieve the same reception pattern by using a suitable output 
filter. ‘This method had the added advantage that the fixed aerial could also 
be used, independently, as an efficient total power system. In this mode of 
operation the spectrum of the aperture distribution is uniform, within a band 
determined by the length of the aerials, and its shape is adjusted to that required 
by a suitable output filter response. The processes of filtering and convolution 
are basically equivalent and in this case where the output of the receiver was 
already recorded in digital form it was convenient to carry out a direct convolution 
with a time function whose Fourier transform provides the required polar diagram. 
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The resulting reception pattern is then identical with that shown in Fig. 2. Noise 
components at frequencies higher than those introduced by the aerial are removed 
before convolution by a simple CR filter. The spectrum of the original signal, 
the convolving function and the CR filter are shown in Fig. 3. 


Fa 


BI 59D 


de. (b) 


Spatial frequency —> 











Fic. 2.—(a) The aperture distribution of the fixed aerial. 
(b) The spatial frequency spectrum of the aerial system. 
(c) The reception pattern. 
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Fic. 3.—The spectrum of the signal (a), the CR filter (b) and the convolving function (c). 


The resultant reception pattern in declination, after the synthesis is completed, 
is that of a 500 ft aperture. The required aperture distribution is here achieved 
by suitably weighting the observations made at the different positions of the 
movable aerial (Ryle and Hewish 1960). 
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(6) The receiver and recorder.—It is necessary, for synthesis, to record quan- 
tities of the form E,E,cos(¢,—¢,) where E,, E,, $, and ¢, represent the 
amplitudes and phases of the voltage induced in the two aerials. This form of 
correlation response can be obtained by the use of a phase-switching receiver 
(Ryle 1952). Inthesystem used here two such receivers were connected to the aerials 
viaa hybrid network as shown in Fig. 4 and arranged to provide outputs in antiphase. 
The difference in output of the two receivers then provides the required quantity 
and the use of the double system improves the signal/noise ratio by V 2 (since 
the combined noise from the two preamplifiers is uncorrelated at the inputs of 
the two receivers). Furthermore, most types of interference are greatly reduced. 


Fixed aerial Movable aerial 
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Fic. 4.—Method of connecting two receivers to the aerials. 


Few sources of interference are received with appreciable intensity by both 
aerials, since the reception patterns do not overlap for a distance of about 20 
miles. Intense and variable interfering signals containing frequency components 
near the switching frequency may however cause a small output even when the 
interference is incident on only one aerial. With the double receiver arrange- 
ment such signals appear in phase at the output and may be rejected with a ratio 
dependent on the gain equality of the receivers. It has been found in practice 
that ignition and lightning interference are attenuated by a factor of 20 compared 
with the output of a single phase- switching receiver which is itself very much less 
susceptible than a simple receiver. 

Since large amounts of data must be recorded | prior to computation, an 
automatic device has been developed for continuously digitizing the signal and 
recording it on punched tape. A normal type of pen recording is also used for 
immediate checking of performance. In the automatic recorder the signal is 
used to operate a self-balancing potentiometer on whose shaft is mounted an 











ior 


fer 
dat 





, i A 


in 
id 


or 








No. 2, 1961 Aperture synthesis method of radio star observation 99 


optical digitizing disk with a cyclic permuting binary code of 8 digits. To mini- 
mize digitizing errors the sensitivity of the receiver is adjusted so that the noise 
level exceeds the digit interval by a factor of 3 or 4. In order to obtain 
synchronism between sampling points made at successive positions of the movable 
aerial the actual sampling rate is chosen to be a sub-multiple of 24 sidereal hours. 
The instant at which the record is sampled (about once every 8 seconds, depending 
on declination) is regulated by a gear train driven from a 50 c/s synchronous 
motor whose power supply is derived from a 100 kc/s sidereal crystal oscillator. 
An impulse from a control circuit disconnects the servo motor in the self-balancing 
potentiometer, thus storing the signal as a shaft position. Subsequent impulses 
cause the outputs of the photo-cells on the digitizer to be read and punched as 
two successive characters using a modified teleprint punch with standard five-hole 
tape. The servo drive is then reconnected and the shaft follows the signal until 
the next sampling time. 

The minimum sampling rate which may be employed is determined by the 
envelope reception pattern (Bracewell 1956) and corresponds to approximately 
four samples of amplitude and phase of the quasi-periodic response in each 
beamwidth (to zero power). In this case the lobes of the interference pattern 
occurred at a frequency slightly greater than twice the minimum sampling rate. 
The output of the receiver was initially sampled four times in each cycle of the 
interference pattern. After further smoothing by the correlation process 
described in Section 2 (a) the sampling interval was increased, during computa- 
tion, to give values of amplitude and phase once per cycle of the interference 
pattern. 

(c) Method of observation.—With the declination of both aerials set to the 
required value the movable aerial is driven to the North end of the rails and 
observations are made for a period of 234 hours. For computation it is convenient 
if all tapes commence at exactly the same sidereal time and an interval must there- 
fore be allowed for moving the aerial and for calibration, etc. before each day’s 
observations. After completion of the first run the aerial is moved a distance of 
42 ft and the observations are repeated for 24 positions of the movable aerial. 
The distance moved is somewhat less than the effective width of the aerial in 
order to eliminate the ‘grating’ side-lobe which would otherwise appear in 
the reception pattern at the extreme edge of the 4°-6 strip (Ryle and Hewish 
1960). It is necessary to maintain an equal path length to the two aerials so 
that a wide bandwidth may be used. ‘This is achieved by a set of cables having 
electrical lengths equal to a whole number of wavelengths and arranged in a scale 
of two; these cables are altered with each aerial move to equalize the path difference 
to within 2-5 wavelengths. 

At low declinations, where the effective North/South aperture is reduced, 
it would be possible to increase the distance between successive positions so long 
as the effective movement in the aperture plane did not exceed about 42 ft. Such 
a procedure would, however, reduce the signal/noise ratio slightly and so far has 
not been employed. 

Throughout the observations separate records are kept of solar activity and 
ionospheric scintillations so that runs may be repeated if serious errors are likely 
to occur in any part of a particular record. Errors due to other sources of inter- 
ference are either immediately evident on the pen record or are indicated by the 
data-checking procedure outlined below. 
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(d) Calibration and checking.—Before proceeding with the analysis of the 
data it is necessary to make a number of preliminary checks. It is particularly 
important to correct for any systematic variations in the sensitivity or phase 
of the system since these will introduce side-lobe responses into the synthesized 
reception pattern in declination. At each of the 24 positions of the movable 
aerial the phases and amplitudes of a number of intense sources are obtained 
from the recorded data. Knowing the accurate declination of these sources it is 
possible to predict the correct value of the phase for each aerial position and so 
to obtain a phase error for each source; provided the sources have a flux density 
greater than about 15 x 10-° w.m~? (c/s)~! they will not be affected appreciably 
by noise or confusion. A mean correction for both phase and amplitude can 
therefore be obtained for each run and used to correct the original data; corrections 
were seldom found to be necessary. 

At the same time as this check is carried out the data are scanned in the com- 
puter for wrong values, caused for example by solar or terrestrial mterference or 
by a fault in the digitizing system. This is done by convolving the data with a 
suitable function and comparing each number on the tape with the ‘expected’ 
value obtained by convolution. Errors detected in this way were reported when- 
ever they exceeded five times the r.m.s. noise level. ‘The necessary corrections 
were then fed into the computer together with the 24 data tapes. The required 
weighting factor for each position, which determines the reception pattern in 
declination, was also fed in at this stage. 

(e) The computation.—It is necessary to carry out the computation in two 
steps since the data tapes must be fed into the computer serially whereas the 
calculation requires information, at a given sidereal time, from all the tapes simul- 
taneously. During the first stage the data are sorted into a convenient layout 
in the store and during this process the convolution mentioned in Section 2 (a) 
is performed. ‘Two convolving functions are used which yield the sine and cosine 
components, sampled once per cycle, of the quasi-periodic output from the 
receiver. ‘These data, corrected for errors as described above, amount to about 
100,000 numbers and are stored on magnetic tape in readiness for the main 
calculation. 

At any given sampling time the data for the 24 aerial positions may be 
represented by the series Cj,, Sjo...... C_12, S_i2. The second stage of 
computation then gives 16 values of D,, 


n= +12 ; 20 

where D,= > (C,+iS,)e" (s = =") > 
n=—12 24 

These values of D represent the output of the full synthesized interferometer, 

sampled at approximately the optimum interval in declination, within a range 

of 4°-6. 

On completion of the calculation it is necessary to correct for the known 
primary reception pattern of the aerials which would otherwise cause the values 
of D at the edges of the declination strip to be systematically smaller than those 
at the centre. While 24 values of C,, and S,, are required only 16 values of D 
are computed since the overlapping of adjacent stations of the movable aerial 
implies that successive values of C,, and S,, are not entirely independent. The 
(complex) values of D are next transformed to |D| and arg D and the printed 
output then gives the 16 values of |D| across the sheet followed by a second line 
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declination running across the page and right ascension downwards. 
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IoI 


|D| is printed as an integer in the range 
o-9999 and arg D as a fraction of a revolution in the range o-g9. The output is 
thus a two-dimensional map of the sky, sampled at discrete intervals, with 


The 


sampling interval is approximately 20’ in declination and 7’-5 in right ascension, 
the precise value of the former depending on the declination of the centre of the 
strip. An example of the printed output is shown in Fig. 5. 
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Fic. 5.—A part of the printed output. 
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Numbers in the upper rows give values of |D| and the 


lower rows the values of arg D; values of D over 99 are printed thus: a = 735. | D|= 71 


corresponds to a flux density cf 107*§ w.m=? (c/s)=}. 
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During the progress of the main calculation a record is made of the probability 
distribution of |D]. This is computed in blocks of 20 successive samples in 
right ascension and the result is obtained on a second output channel. This 
information is required for statistical investigations of the distribution of radio stars 
(Scheuer 1957) and the application of the method will be discussed in another paper. 

(f) Analysis and interpretation of the results.—The position of a given source 
in right ascension is determined precisely by the value of the phase (arg D) 
but it is necessary to interpolate, both in declination and right ascension, to find 
the accurate declination and flux density of asource. Since the reception pattern 
is known in both coordinates this interpolation is most conveniently carried out 
by a simple graphical process. In cases where the envelope reception pattern 
appears to have the wrong shape, it is probable that more than one source is 
present and values of the phase then give valuable information which may resolve 
the blend. 

No attempt is made to reduce the collimation errors to zero for each series of 
observations. Instead, use is made of the accurate positions of 64 intense radio 
stars which have already been determined in an independent investigation 
(Elsmore, Ryle and Leslie 1959). Any declination strip will usually contain 
several of these sources and it is therefore possible to obtain both collimation 
errors and a flux density normalizing factor for each series of observations. 


3. Performance. 


(a) The reception pattern.—It is only necessary to consider the reception pattern 
along lines in right ascension and declination which intersect at the centre of the 
pattern. At other angles the response is given by the product of these patterns 
and side-lobes will be very small. 

The observed reception pattern in right ascension, derived from a uniform 
excitation of the fixed aerial and the convolution method described in Section 
2 (a), is shown in Fig. 6 (a); it is seen that the side-lobes are less than 0-5 per cent 
at angles greater than 1°. 
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Fic. 6.—The observed reception pattern (a) in right ascension and (b) in declination. 
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In declination the reception pattern will be affected by errors which may occur 
in the phase or amplitude of the signal. While systematic errors can be removed 
as discussed in Section 2(d) there will be random errors due, for example, to 
refraction inthe ionosphere. Variations of amplitude are found to be of negligible 
importance compared with the phase variations which have an observed r.m.s. 
value of about 15°, corresponding to a figuring error of about A/25 in the synthesized 
aperture. The latter introduce side-lobes with an r.m.s. level of about 6 per cent 
confined to the 4°-6 envelope; beyond this the response is less than 0-5 per cent. 
The observed reception pattern in declination is shown in Fig. 6 (6). 

The cause of the random phase variations is not yet known. From existing 
knowledge of non-uniformity in the ionosphere it seems unlikely that they can 
be explained simply by ionospheric refraction and it is possible that they may 
be due in part to temperature variations in coaxial cables. 

(b) Sensitivity.—The noise level for the initial observations at declinations 
of 50° and 52° was about three times higher than that of the later observations 
owing to errors introduced by the digitizer. For the surveys at declinations 
of 42° and 05° an improved digitizer was used and uniform excitation of the long 
aerial was adopted. ‘To distinguish between true noise and the effect of weak 
unresolved sources, special records were made with one of the aerials replaced 
by a resistance at a temperature close to the effective sky temperature. Under 
these conditions there is no response from radio stars but the noise level due to 
the combined effects of preamplifier noise and the galactic background emission 
is the same; additional noise components due to weak terrestrial interference 
incident on one aerial are also included. The possible incidence of coherent 
interference in both aerials was investigated in a separate experiment and shown 
to be unimportant. 
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Fic. 7.—The integrated probability distribution of |D| for the 42° survey —+—*—, and that 
derived —o—0o— from the noise measurement. 
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The observed noise level was found to be 0-08 x 10-*6 w.m~? (c/s)—! as com- 
pared with an estimated theoretical value of 0-045 x 10-® w.m-? (c/s)-1. 

The integrated probability distribution of |D| for the 42° survey is shown in 
Fig. 7 together with that derived from the independent noise measurement. It 
is evident from these curves that the noise on the records is likely to be unimportant, 
either in the measurement of individual sources, or in the statistical analysis of 
|D. 


4. Observational results 


Four surveys have been made centred on declinations of 52°, 50°, 42° and 
05°. The first two (52°, 50°) were carried out primarily to test the operation 
of the whole system. A declination strip 2°-5 wide was included in both surveys 
and a direct comparison of the results has been made. The noise and side-lobe 
level were both higher than in the later surveys which used an improved digitizer 
described in Section 2 (6). 

The survey at 05° was made primarily to investigate the distribution of sources 
for different directions in the galaxy. The most recent survey at 42° represents 
the best series of observations for the study of individual sources and for application 
of the statistical method of analysis. 

(a) Comparison of surveys at declinations of 52° and 50°.—The results obtained 
for the two surveys are entirely independent, both with respect to errors in flux 
density and position and to those arising from side-lobes. A comparison of 
results in the overlap region 2°-5 wide in declination thus provides an excellent 
overall test of the operation of the instrument. 

In Table [ a list is given of all sources in the overlap region having flux densities 
greater than 3 x 10-*° w.m~* (c/s)! and it is seen that there is extremely good 
agreement. Of the 67 sources listed, the majority show good positional agree- 
ment and have flux densities which agree to within 20 per cent. ‘Two sources 
found on the 50° survey had flux densities of less than 2 x 10-76 w.m-? (c/s)~! on 
the 52° survey; in one case the wrong right ascension lobe was chosen. 





TaBLe I 
52° Survey 50° Survey 

R.A. Dec. Flux density R.A. Dec. Flux density 
h m $ ° 10-*%* MKS h m - . 10-28 MKS 
eo 13 213 50 |22°4 5°7 00 13° 21°0 5° 192 4°0 
Co 35 43°5 50 617°6 50 Co 35 642°9 50 6166 5°6 
©0 40 214 51 46°7 37 co 40 19°9 51 47°5 35 

<7'5 00 47 #475 51 42°3 471 
co 48 06:2 50 656-7 11°4 co 48 05°6 50 561 Ir‘! 
Ol 19 44°77 49 51°7 4°3 Or 19 455 49 49°8 41 
02 23 38:6 51 0o6°5 3°5 02 23 «381 51 03°6 3°8 
02 40 175 49 51°4 5"0 02 40 190 49 538 5°5 
03 «28 =«20°1 51 25°9 3°3 03 28 «19°7 Sz 2a°7 3°1 
03 34 «=(O4'! so 37 12°0 03° «34 «=(OF4'l 5° 35°2 12°5 
03 «55 47°6 50 6488 6-1 °3 55 46:0 5° 49°7 6°5 
04 26 o1°2 5° 093 3°9 04 25 581 50 I0°1 43 
04 37 369 51 25°7 24 04 37 31°2 51 25°0 3°2 
04 42 564 50 32°0 69 04 42 561 5° 33°5 73 
04 48 55°6 52 000 6°5 04 48 533 52 o2'1 5‘1 

<2°0 04 56 185 50 092 $3 
©5 12 41°6 51 04% 4°8 ©5 12 39°0 51 097 4°8 
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™ on? Qaseiny TABLE I—continued of? Cee 
R.A. Dec. Flux density R.A. Dec. Flux density 
n h m 8 4 ? 10-** MKS h m 8 ° ’ 10-7 MKS 
it 05 15 40°2 50 49°8 9°0 05 15 381 50 52°6 8-9 
t 05 22 I1'0 51 37°4 6-2 05 22 085 51 39°1 5°9 
f 05 38 446 49 487 53 05 38 42°9 49 510 64 
06 «15 587 51 03°2 2°9 06 #15 57°0 51 088 2°7 
06 25 54°0 50 628-6 4°3 06 25 53°9 50 32°2 4°9 
06 28 49°9 50 23°3 3°1 06 28 = 50°2 50 27°4 2°5 
06 35 «45°2 51 085 3°2 06 35 44°8 51 084 2'1 
06 50 15°2 50 «161 4'1 06 «50 ~= «6186 50 24°38 2'0 
id 07 00 30°! 50 3370 371 07 00 31°6 50 35°7 2°5 
mn 07 03 47°99 50 04°5 4°9 07 03 465 50 0973 4°3 
ys oT a 17% 5° 393 3°0 07 24 16°6 50 42:2 2°9 
ye 07 28 37:9 51 17°6 3°1 o7 28 357 5I 20°0 2°6 
“ 07 29 6487 51 58-2 4°2 07 29 45°5 52 O1'9 2°7 
07 36 30°1 51 30°0 4°6 07 36 27°9 51 39°0 4°6 
08 29 44°6 51 096 6:9 08 29 404 5I 173 5°8 
es 08 35 33°8 51 096 3°5 08 35 31°0 51 16:0 2°1 
ts 08 35 587 50 29°6 3°4 08 35 55°2 5° 35°9 2°5 
n 0&8 39 394 51 O4'l 3°2 08 39 344 SI 114 2°4 
09 00 25°4 50 24'1 3°1 ©9 0O 21°7 50 30°1 3°2 
10 31 208 5° 29°5 6:0 sO 3% 17% §0. -31°3 5°9 
ed 10 39 089 =: 30°9 5°6 10 39 077 50 293 5°5 
ux Il 30 347 50 2773 3°2 Ir 30 313 50 2971 41 
of tr 36 13°0 50 36-2 3°2 1m 96 «ra 50 35°7 2°6 
nt Ir 43 078 50 044 9°4 II 43 05°9 5° 023 8-9 
II 50 52°09 49 5§2°3 5" Ir 50 498 49 50°8 5°5 
: oe ae eg 31 II 51 105 $1 34°2 3°4 
es 12 20 296 50 49°1 3°7 I2 20 27°0 50 47°4 3°6 
od I2 49 29°09 50 52°4 4°3 12 49 285 50 52°1 4°5 
“e- 13 17 44'5 52 06:2 2°9 13 17 41°7 52 063 3°5 
“es 13 290 371 50 261 41 13 29 39°90 50 25°! 4°9 
on 13 43 202 50 O35 | 7% 13 43 292 50 009 11'S 
14 41 244 50 5I°7 4°6 14 41 21°1 5° 49°9 3°3 
14 46 59°5 5° 30°5 2'8 14 46 59°7 5° 25°5 3°2 
14 52 24°7 50 12°7 4°2 14 52 248 50 I5°0 4°7 
16 52 40°5 52 05°2 3°5 16 52 43°6 52 03°2 3°7 
17 00 102 50 49°5 3°2 17 00 13°32 50 543 3°2 
sity 17 23 06°0 51 00°5 I1°2 17 23 «21178 5I 02% 10°8 
KS 17 29 51°6 5° 05°7 3°8 17 29 56°7 50 I1°4 471 
16. 17 -. So's 50 28'5 74 18 17 568 50 6281 7°6 
19 40 21°4 50 «281 8°5 19 40 26°7 5° 29°5 7°93 
20 OF 49°4 50 08-6 5°3 20 O07 54°4 50 12°1 38 
20 30 12°5 51 45°6 5°9 20 30 17% 51 45°3 5°7 
20 37 O7°1 51 o7'1 12°5 20 37 104 5I ogo I2‘I 
20 40 363 50 17°2 6°7 20 40 42°2 50 21'1 71 
20 55 446 50 47°1 4°3 20 55 49°5 50 481 4°5 
21 22 43°0 5° 49°4 3°5 21 22 46:9 5° 45°4 3°7 . 
at 23 639 51 54°2 4°2 21 23 006 St ‘§7°t 4°5 
21 50 15°4 49 45°3 4°2 Region not covered by survey 
23 45 343 51 056 3°3 23 45 362 51 03°9 3°0 
23 46 088 50 26-4 3°3 23 46 55°09 50 17°9 3°0 
A precise comparison of the positions yields an r.m.s. differential error of 
~ 3’ in declination and ~ 48in right ascension. The declination errors appear to be 
random but there is a systematic error in right ascension which, when plotted 
against right ascension, approximates to a sinusoidal variation as shown in Fig. 8. 
8 
; 
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Fic. 8.—A plot showing the systematic error in right ascension between the 50° and 52° surveys. 


The r.m.s. error of the individual surveys thus amounts to ~1’ in right 
ascension and ~2’ in declination. When the systematic error is more fully 
understood the error in right ascension should be reduced to ~0’:3. 

(b) Survey at o5°.—Previous observations have indicated that while most 
radio stars remote from the galactic plane are probably of extra-galactic origin 
there is a possibility that they can be accounted for by a class of objects in the 
galactic halo (Ryle 1958). This possibility may be eliminated by studying the 
isotropy of the distribution of sufficiently weak radio stars. A survey has there- 
fore been made at a declination of 05° which includes longitudes towards the 
centre and anti-centre at galactic latitudes near 50° where the greatest anisotropy 
might be expected to occur. 

Application of the statistical method enables deductions to be made about the 
number-intensity relationship of radio stars at a flux density considerably lower 
than that at which the sources can be reliably distinguished individually. The 
probability distribution of |D| was computed for two areas towards the centre 
and anti-centre (right ascensions 22? 40™—o02"40™ and 10% 40™— 14 4o0™), 
The results are shown in Fig. g. It is seen that any differences in these distri- 
butions are less than the statistical uncertainty and hence there is no evidence 
for a reduction in the density of radio sources in the longitude of the anti-centre. 
It will be shown in another paper (Hewish 1961) that this observation implies 
an isotropy which extends to at least 10,000 sources steradian. 

(c) The 42° survey.—Positions and intensities were derived for sources having 
flux densities >2 x 10-*6w.m~*(c/s)-4. 175 sources were found in an area 
of 0-3 steradians («= 20" 40™—19%15™, 5=40° — 44°); this number corresponds 
to an average of 1 source per 30 beam areas so that errors due to confusion should 
be small. A complete list of these sources is given in Table II. The quoted 
error for each source includes a confusion error derived from the probability 
distribution of |D|, a reading error and the residual collimation error; it does not 
include the systematic error in right ascension, discussed in Section 4 (a), for 
the 50°— 52° strip, which has an upper limit of + 4°. 
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The radio sources M31 and NGC 1275 came within the region of the survey. 
A number of weak sources were found within the boundaries of M31 but, as 
discussed elsewhere (Baldwin and Costain 1960), there is no significant excess as 
compared with a neighbouring area of sky and hence no reason to associate them 
with the nebula. NGC 1275, a member of the Perseus cluster, has been resolved 
into two distinct components separated by 32’, together with a weak extended 
component which may be associated with the cluster as a whole. The new 
information about NGC 1275 will be discussed fully in another paper (Leslie 
and Elsmore 1961). 
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Fic. 9.—The probability distribution of |D| for areas at galactic latitudes near 50° and longitudes 
towards the centre and the anti-centre. 
@ R.A. 222 40M to o2h4om 
x R.A. 102 40M to 14h 4om 


Tasie II 

R.A. Error Dec. Error Fluxdensity Error 
ee ig +s 2 , +° 10" MKS) s+ 
00 IO 54°7 3°0 40 34°6 3°0 5°0 06 
00 19 09°6 2'0 43 10°7 3°0 3°5 06 
00 20 50°0 2°5 43 48-0 3°0 2°5 o'5 
00 28 06°7 a°s 40 53°2 3°0 2°5 o's 
00 41 54°I 2°5 42 30°5 3°0 2°0 04 
00 45 25°6 2°0 40 06:0 3°0 3°2 05 Possibly complex 
00 51 38-4 2'0 40 26°1 2°5 3°2 o's 
00 53 581 2°5 43 58:0 830 2°4 0'5 
OI 03 51°6 2°0 42 128 2'0 K te J o'5 
OI og 22°7 20 41 37°0 2°5 3°6 o'5 
Or 13 33°8 30 409 03°3 3°0 2°3 0°4 Possibly complex 
OI 20 32°0 a°s 40 34°8 4°0 28 o'5 Complex 
Or 44 5532'S 43 09°3  _3°0 2°5 0°5 
Ol 52 263 Is 43 32°5 2°0 9°6 1°0 
ol 53 17°3 2°0 41 48-9 2°5 31 0°5 
Or 57 03°3 2°5 49 34°! 3°0 2°7 o's 
02 16 00°6 2'0 42 24°4 2°5 3°9 06 
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h m 8 
02 20 083 


02 22 37°0 
02 32 45°6 
02 46 15°0 
02 54 33°9 
03 00 16°! 
03 02 30°! 
03 «09 «42°6 
03 +%&II 23°6 
03 14 54°2 
03 16 28°6 
03 «23 «2183 
o3 37 607° 
03. ge §°°7 
03 44 37°8 


03 54 19°6 
04 04 35°7 
04 16 49°! 
04 19 52°6 
04 22 37°8 
04 25 47°2 
04 29 08:0 
04 41 #17°7 
04 46 04°4 
O§ Ir 12°! 
05 25 15°5 
05 28 I0°1 
06 02 43°0 
06 II 24°0 
06 12 §2°5 
06 13 32°5 
06 14 47°3 
06 21 36°5 
06 21 50°! 
06 27 13°3 
o6 «28 46°71 
06 42 51°8 
06 44 12°! 
06 47 22°4 
06 49 29°! 
06 52 37°4 
07 03 104 
07 09 02°6 
07 26 19°0 
oF 3 33°7 
08 of 42°7 
08 06 39°5 
08 09 33°4 
08 20 O5°1 
08 36 39°0 
08 36 56°0 
08 49 17°7 
08 56 47°5 


09 OI 002 


Error 
+s 
1°5 
2'0 
1°8 
1°5 
2°5 
2°5 
2°5 
2°s 
1°5 
2'0 
1°5 
1°8 
1°8 
2°5 
3°0 
2°5 
2°5 
1°5 
2°5 
1°8 
2°5 
1°8 
1°5 
2°5 
2°5 
2°0 
1°8 
2°5 
20 
2'0 
2°5 
2°5 
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TaBLe [I1—continued 


Dec. 


49°6 
19‘2 
II‘l 
56-2 
40°4 
35°5 
55°5 
24°7 
06°0 
44°7 
21°0 
44°6 
53°4 
24°4 
39°1 
47°2 
42°4 
54°1 
37°2 
47°7 
21°2 
27°7 
33°7 
54°0 
13°0 
37°5 
40°5 
46°1 
37°9 
53°1 
31°4 
41°4 
33°5 
O1'3 
30°6 
31°4 
20°3 
49°7 
22'9 
33°5 
39°6 
42°2 
35°3 
55°3 
08-0 
49°2 
44°5 
36°8 
26°4 
07°6 
37°9 
13'0 
30°2 
34°6 
50°0 


Error 


+’ 
2°0 
2°5 
2°0 
1°8 
3°0 
3°0 
3°0 
3°0 
1°8 
3°0 
1°5 
2'0 
2'0 
3°0 
3°0 
3°0 
3°0 
I°5 
3°0 
20 
3°0 
2°0 
1°5 
3°0 
3°0 
2°s 
2'O 
3°0 
3°0 
2°5 


Flux density Error 


10-76 MKS 
16°6 
3°1 
4°5 
6°6 
2'0 
2'1 
2°8 
21 
8-5 
14°5 
41°0 
5°2 
5°5 
3°0 
2°8 
2°5 
2°4 
25°1 
2°8 
51 
2°8 
4°2 
14°9 
2°I 
2°I 


x 
1°5 
o'5 
o's 
06 
04 
o*4 
o'5 
o*4 
o°8 
2°0 
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Possibly complex 


Complex 
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TasLe II—continued 
R.A. Error Dec. Error Flux density Error 

. “a. e +s ° ’ =" 10° MKS + 
09 02 50°0 2°5 41 30°2 3°0 3°0 06 
09 03 50°2 4°0 42 55°4 5"0 2°7 0°5 
09 04 22°2 2°O 41 49°4 2°5 5°2 "7 
09 06 19°5 1°5 43 05°0 2°0 19°5 1°5 }Confused region 
09 06 34°8 5°0 40 46:2 5'0 4°5 1‘o 
09 06 46°6 5°0 41 190°7 5°0 4°0 o'8 
og 22 I4'! 3°5 42 31°6 5'0 2°4 0'5 
09 aa §2°2 3°5 4a 27'3 50 2'0 o'5 
09 35 06°7 2°5 42 53°38 3°0 2°5 0"4 
09 36 14°0 2"0 40 30°2 2°0 3°4 0°5 
09 45 25°2 1°8 41 53°6 2°0 4°I 06 
09 45 53°0 2°5 40 47°6 3°0 2°3 o"4 
09 47 «+I3°! 2°5 42 28°0 3°0 2°0 0°4 
10 07 27°4 I°5 41 49°9 1°8 6-0 o'8 

10 09 08°7 °5 43 30°8 3°0 2‘1 o'4 

ex 
10 13 008 2'0 40 59°6 2°5 3°8 06 
10 20 14°6 2°5 40 05°4 3°0 2°3 0"4 
10 22 33°0 2°0 43 15'2 2°5 4°5 o'7 
10 33 27°0 3°0 40 48°7 5°0 2°0 0-4 
10 34 20°4 2°5 49 27°7 3°5 371 o°5 
10 55 10°2 3°0 40 27°8 4°0 2°3 o"4 
10 56 10°0 1°5 43 17°8 Ts 11‘2 I°4 
11 08 53°! 2'0 41 08°4 2°0 5'1 o'7 
II 09 53°4 2'0 43 43°! 2'0 63 08 
II II 54°0 I°5 40 52°0 1°5 21°2 I's 
11 28 07°4 2°5 43 42'2 4°0 3°2 06 
it 3% 58-5 1°5 43 47°4 1°8 6°6 o'8 
II 32 06°2 3°0 40 56°6 3°5 2°8 o's 
12 06 45°0 1°5 43 56°6 1°5 10°7 1°0 
12 20 08:8 2°5 40 50°9 3°5 2°O 0°4 
12 22 03°0 I°5 42 22°7 18 7°2 o°8 
12 31 5§5°0 2°5 43 10°7 3°0 2°5 0°5 
3 32 of6 2°5 41 26°5 3°0 2°4 o*4 
12 49 44°7 2°5 43 18-6 3°0 2°3 0"4 
3 83 169 2°0 43 13°6 2°5 a°s o'5 
12 58 14°6 1°5 40 23°0 1°8 8-6 1‘O 
13 19 O81 2'0 42 51°6 2°0 6-0 08 
13 19 41°7 4°0 43 27°7 5°0 2°4 0°5 
i 2t 14's 30 41 30°0 3°5 2°9 o'5 
13 24 54°! 2°5 43 O81 3°0 2°1 "4 
13 33 «+10°7 1°8 41 148 2'0 5‘0 o'7 
13, 40 27°3 3°0 43 51°5 3°0 2°2 o'4 
13 43 «28°97 2°0 43 04°6 2°0 3°4 06 
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5. Conclusion.—The first surveys carried out at a wavelength of 1-7 m using 
the new principle of aperture synthesis have shown that the theoretical possibilities 
of the technique may be fully realized in practice. 

A primary resolving power of 25’ x 35’ has been achieved with a noise level 
of 0-08 x 10-*6 w.m~*(c/s)~-!. The side-lobes in right ascension are negligible 
beyond +1° from the centre of the beam; in declination the side-lobe level is 
~5 per cent within +2}° and less than o-5 per cent at greater angles. Radio 
stars may be reliably detected to a flux level of 2-0 x 10-* w.m~* (c/s)~, —1, corres- 
ponding to one source/20 beam areas, with positional accuracies of +2’ in 
declination and +1’ in right ascension; at this limit the noise is still only 5 per 
cent so that statistical investigations may be used to study the distribution of 
radio stars at smaller flux densities. 
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OBSERVATIONS OF THE ZODIACAL LIGHT FROM A VERY 
HIGH ALTITUDE STATION 


I. THE AVERAGE ZODIACAL LIGHT 
D. E. Blackwell and M. F. Ingham 


(Received 1960 October 6) 


Summary 


The paper presents measurements of the absolute surface brightness at 
three wavelengths, the polarization at two wavelengths and the position of 
the zodiacal light made from the high altitude station of Chacaltaya in 
the Bolivian Andes (height 17 100 ft, geomagnetic latitude 3°S). This site 
is particularly suited to such observations because the weather is good and 
the atmospheric extinction is close to that expected for a dust-free Rayleigh 
scattering atmosphere. The brightness and polarization measures extend 
to within 19° of the Sun. The colour is close to that of the Sun, but it is 
slightly more red; the polarization measures agree well with those recently 
obtained by Elsasser. 





Introduction.—The zodiacal light has been the subject of many investigations 
during the last three centuries, although accurate measurement using modern 
techniques may be said to have begun only twenty years ago with the work of 
Elvey and Roach (1). It is now recognised, as a result of these investigations, that 
most of the zodiacal light originates in interplanetary space, and that it may be 
considered as the outermost part of the solar corona. Also, it has been supposed 
that much of this extra-terrestrial part, and perhaps all, is due to the scattering of 
sunlight by small solid particles, usually described as ‘‘dust’’, situated in inter- 
planetary space at distances greater than 0-3.A.U. from the Sun—as was first 
suggested by Cassiniin 1672 (2). This viewpoint was almost universally accepted 
until 1953, when Behr and Siedentopf (3) made the first accurate measurements of 
the degree of polarization of the zodiacal light, and on the basis of these measures 
put forward the hypothesis that approximately half of the surface brightness of 
the zodiacal light near elongation 35° is due to the scattering of sunlight by free 
electrons in interplanetary space, the electron density in the plane of the ecliptic 
at 1A.U. from the Sun being about 600 electrons cm-*. ie remainder of 
the zodiacal light is supposed due to scattering by interplanetary dust particles. 

It has seemed to the authors that there are good arguments against the existence 
of such a high concentration of free electrons in the plane of the ecliptic so far from 
the Sun (see, for example, Blackwell (4) and later papers in this series) and that the 
presence of these electrons should be confirmed by another and independent 
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method. This is the origin of the present investigation, which was made during 
an expedition planned for this purpose and sent from the Cambridge Observatories 
to the high altitude mountain site of Chacaltaya in Bolivia during 1958. As part 
of the programme of the expedition we took advantage of the exceptionally fine 
atmospheric conditions at Chacaltaya to re-evaluate the fundamental data con- 
cerning the brightness, polarization and position of the zodiacal light, and studies 
of possible variations of these quantities ; these studies are described in two papers 
of this series. In a third paper we shall describe the spectrum of the zodiacal 
light that was obtained at Chacaltaya and discuss its bearing on the problem of the 
existence of interplanetary electrons, usually known as interplanetary gas. Ina 
fourth paper we shall rediscuss the properties of the interplanetary dust. 

Observing site.—The Chacaltaya observations described here have been made 
only on the inner zodiacal light within the range of elongation, 19°<«<70°, 
There are two reasons for this restriction ; the first is that the problem of detecting 
free electrons in interplanetary space becomes rapidly easjer as the elongation of 
the line of sight diminishes—for instance, at elongation 19° the line of sight 
approaches to within 0°33 A.U. (= 70R,) of the Sun—and indeed, the observing 
site was especially chosen so that observations could be made as close to the Sun 
as possible. The second reason is that at greater elongations there is a large and 
very uncertain correction for the night sky background, which requires special 
methods for its evaluation, and which in any case cannot be better investigated at 
Chacaltaya than at any other suitable mountain site. 

The requirements for a site suitable for observations such as these are; (a) 
a situation in the tropics, from which the zodiacal light may be observed when its 
symmetry axis is nearly vertical at the horizon, (6) a very transparent atmosphere, 
so that photometric observations may be made with reasonable accuracy (i.e. with 
an error not greater than 5 per cent) to within a few degrees of the horizon, and (c) 
a good climate which enables observations to be made near to the horizon on 
successive nights with little interruption by cloud. The cosmic ray station of 
Chacaltaya, directed by Professor I. Escobar V, was found to fulfil these require- 
ments well. The station is in the Bolivian Andes at a height of 17 100ft, at 
latitude — 16° 19’, longitude 68° 10’, and geomagnetic latitude — 3°. It possesses 
a good west horizon at zenith distance 91° where the sky was usually of extreme 
clarity at sunset and during twilight and night, although the weather during the 
day was frequently poor. Our observing programme was scarcely interrupted by 
cloud. . 

As part of the purpose of the expedition was to obtain the absolute brightness 
of the zodiacal light and to seek a variation of brightness with time, the atmospheric 
extinction was measured whenever the zodiacal light was observed. The 
extinction measurements were made photoelectrically at the two effective wave- 
lengths of 0-45 4 and 0-59, using bright stars in the region of the zodiacal light 
and a distant standard lamp to keep a check on possible variations in the sensitivity 
of the photometer. A logarithmic plot of one set of observations at 0-45 1, which 
is typical of all, is given in Fig. 1; here the brightness is related to the air mass 
calculated for a spherical earth, and the regularity of the observations show that 
the atmosphere is reasonably uniform within to a few degrees of the horizon. To 
illustrate the quality of the atmosphere at this site we give in Table I the atmo- 
spheric transmission measured on 23 evenings at zenith distance 75° over the 
western horizon (corresponding to an air mass of 2-35) for wavelengths 4500A 
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and 5900A. ‘The zenith absorption at 5900 A is usually about o™-08. ‘The data 
are summarized in Table II, from which it will be seen that the extinction at both 
wavelengths is close to that expected theoretically for a dust-free Rayleigh 
scattering atmosphere (5). 
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Fic. 1.—Variation of the brightness of a star with path length. 


TABLE | 


Atmospheric transmission (per cent) at zenith distance 75° at Chacaltaya (altitude 17 100 ft.) 


Date 4500A 5900A Date 4500A 5900A 
June 4 60°0 72°0 July 2 60°0 72'0 
5 58-0 69°6 3 61°0 73°0 
9 55°8 66°8 4 58°9 70°7 
10 59°0 72°0 5 51°0 61°1 
12 54°5 65°3 7 520 62°0 
13 52°9 66-9 8 50°3 60°6 
15 60°0 69°8 9 49°2 59°0 
16 57'°5 65°4 II 600 72'0 
17 60°0 73°2 14 56-7 68-0 
18 60-0 72°0 15 5673 66°5 
19 55°9 69°8 16 57°0 68°4 
17 61°0 73°0 


Apparatus.—The surface brightness and the degree of polarization of the 
zodiacal light were measured photographically using three cameras. A central 
camera, which was also used to measure the position of the axis of symmetry of the 
zodiacal light, was pointed directly at the ecliptic, and two others photographed 
areas of the sky near the horizon on each side of the ecliptic. All the cameras had 
simple lenses of 7 cm diameter and relative aperture F/1-7, giving a plate scale of 
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TaBLe II 
Summary of atmospheric transmission data 


Percentage transmission 


4500A 5900A 
Transmission at z=75° (average) 54°2 68-1 
(maximum) 610 73°2 
Zenith transmission (average) 86-1 QI'5 
corresponding absorption (o™-162) (0™-097) 
(maximum) 88-9 93°I 
corresponding absorption (o™-128) (o™-078) 
Exponential absorption coefficient corresponding to maximum 
transmission o°210 0°133 
Theoretical absorption coefficient for Rayleigh scattering (5) o'219 0°134 


1°=2:14mm._ The central camera used four lenses and photographic plates of 
size 12in. x 12in. on which the four images were simultaneously exposed, the field 
covered by each individual camera being about 56° in azimuth and 50° in altitude. 
The side cameras, containing two lenses each, photographed somewhat smaller 
fields. ‘Two of the lenses of the central camera were used to photograph the 
zodiacal light simultaneously through the Chance filters OBr1o (blue) and OY1 
(orange) on Kodak IIaE plates, and also sometimes through the Chance filter 
OV1 (ultra-violet) on IlaO plates. The other two lenses of the central camera 
were used for the simultaneous measurement of polarization by a method which 
will be described later. The side cameras were used with the same filters as the 
central cameras for the measurement of brightness only. As the exposures were 
usually of about 10 minutes’ duration, no guiding was necessary. After the 
photographs had been taken, calibration marks were impressed upon previously 
shielded parts of the plates using a step wedge and standard lamp, reduced in 
brightness with suitable neutral filters, and all the plates, from the central and 
side cameras, were processed together. 


Measurement of the brightness of the average zodiacal light.—In this paper we 
take no account of possible variations of the zodiacal light, but merely average the 
data to give an observational model representing the ‘‘average’’ zodiacal light. 

The surface brightness of the zodiacal light has been measured before on 
many occasions but, so far as we are aware, never under such good atmospheric 
conditions as exist at Chacaltaya. One of the difficulties in work of this kind is to 
distinguish between the zodiacal light and the background light of the night sky, 
and imperfect separation of these two is probably the chief source of uncertainty 
in the best of the modern data. A particular advantage of attempting this separa- 
tion at a site such as Chacaltaya is that there the observed ratio of zodiacal light to 
night sky brightness is rather greater than usual because the atmospheric extinction 
is small and because observations can be made at small elongations where the 
zodiacal light is bright. We use the term ‘‘night sky”’ in a rather loose sense 


here; the zodiacal light could usually be seen when the Sun was depressed only 
12° below the horizon and some of our observations of brightness (although none 
of the spectra) were made before the end of astronomical twilight. 

A typical photograph taken with the central camera is shown in Plate 3; this 
photograph was taken through the OB10 filter using an exposure time of 10 minutes, 
the Sun being 19° below the horizon at the mid-point of the exposure. The 
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quality of the atmosphere is demonstrated on the photograph by the presence of 
star trails which may be seen intersecting the depressed horizon at a distance 
of 80 miles. 

We anticipate the detailed results here to give the response of the cameras to 
the zodiacal light as a function of wavelength. By this term we mean the product 
of the following factors; transmission of filters, transmission of camera lens, 
sensitivity of photographic emulsion, emission of zodiacal light, atmospheric 
transmission. ‘The transmission of the camera lenses, which were made of white 
glass, was determined with a Carey spectrophotometer, and to make this easier 
an optical flat was polished on each side of alens. ‘The transmission of each filter 
was measured to an accuracy of 1 per cent of the incident light. ‘The response 
curves for the three filters are shown in Fig. 2. 
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Fic. 2.—Response curves for the three filters (the response is the product of filter transmission, 
transmission of camera lens, sensitivity of photographic emulsion, emission of zodiacal light, 
transmission of Earth’s atmosphere). 


The effective wavelengths at which these three filters measure the emission of 
the zodiacal light depend upon the energy distribution in the spectrum of the 
standard lamp, and we simply refer here to a description of a method of deter- 
mining these wavelengths (6). The three wavelengths are given in Table III. 


TaBLe III 
Filter Effective 
wavelength 
OV1 3880A 
OB1o 4470 
OY1 6200 


The standard lamp was a 4v. 4 watt lamp selected for the optical quality of its 
envelope and run at the colour temperature of 2285 °K. It was calibrated by the 
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National Physical Laboratory, Teddington, over the range 4000 A to 7000 A with 
an error of not more than 1-2 per cent, and with rather more uncertainty at the 
shorter wavelengths. 

Reduction of photographs.—To obtain the combined brightness of the zodiacal 
light and night sky background, the plates were scanned parallel to the horizon 
(i.e. approximately perpendicular to the ecliptic), using the Cambridge Observa- 
tories’ recording microphotometer with a scanning spot of size 0-5 mm x 0°5 mm. 
The photographic densities so obtained were converted to brightness values using 
the calibration, and these energies were then multiplied by a camera vignetting 
factor which had been measured photoelectrically in the laboratory. A typical 
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Fic. 3.—Typical scan across the zodiacal light reduced to intensities. 


scan across the zodiacal light, at elongation 45°-5 and altitude 24°-3, which shows 
the brightness out to 40° on each side of the ecliptic is reproduced in Fig. 3. Such 
a scan results from combining microphotometer scans of the central photograph 
with scans of the two side photographs. 

From scans such as these the ratio of the zodiacal light brightness on the 
ecliptic to the background brightness has been deduced using the methods and 
results of Roach (7). At these smaller elongations the error due to uncertainty 
in the brightness of the sky background is not great, being no more than 5 per cent. 
The observed zodiacal light brightness is now corrected for atmospheric extinc- 
tion taking the altitude of the line of scan from the photograph. However, the 
correction for extinction is not straightforward. Blackwell (4) has suggested that 
it is not sufficient to allow for extinction by a simple application of the usual 
exponential relation with an absorption coefficient determined from observations 
of stars. ‘The reason for this is that the zodiacal light is an extended source, and 
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some of the light from the outer parts will be scattered by the Earth’s atmosphere 
into any point of observation, producing a spurious component. ‘The effect has 
been examined quantitatively by Fessenkov (8) taking an extinction law of the form 


a sec z 
Loos 7” Tinciaent P 


and putting p equal to 0-83, corresponding to an exponential coefficient of 0-184. 
Fesenkov estimates the amount of spurious light to be 5-8 per cent on the axis of 
the zodiacal light at zenith distance 75°. At Chacaltaya, because of the lower 
atmospheric pressure, the amount of spurious light on axis may be expected to be 
approximately half of this (3 per cent) and to be less still in the red region (2 per 
cent). We do not regard the correction as an important one in measurements of 
absolute brightness at smaller elongations along the ecliptic, and especially if the 
observations are made at a very high altitude station, although, as Fesenkov 
points out, it is an important correction to measurements away from the axis. 

In relating our absolute measurements of brightness to the average brightness 
of the solar disk, we use the data of Johnson (9g) for the solar energy distribution. 

Results of measurements of brightness.—The results of these measurements are 
shown in Figs. 4 and 5 where the brightness of the zodiacal light along the ecliptic, 
corrected for extinction, is plotted as a function of elongation for the two principal 
wavelengths used; the unit of brightness is that of the integrated solar disk. 
These diagrams show the brightness data for 22 evenings between the dates 1958 
June 5 and Aug. 2, each night being represented by at least four points, and more if 
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. 4.—Brightness of zodiacal light at wavelength 6200A, OY 1 filter 
(unit; brightness of the integrated solar disk). 
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one photograph is available. 
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Fic. 5.—Brightness of the zodiacal light at wavelength 4470A, OB1o filter 
(unit; brightness of the integrated solar disk). 
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The spread of the results is rather great, 
being about 15 per cent at elongation near 35°, and we shall consider in Paper II 
of this series how far this spread can be interpreted as a real variation in brightness, 
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Fic. 6.—Distribution of energy in the spectrum of the zodiacal light in comparison with that 
of the integrated solar disk. 
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Fig. 6, obtained from the present observations, shows the distribution of energy 
in the spectrum of the zodiacal light at elongation 35°, corrected for extinction, in 
comparison with that for the average solar disk ; this diagram also includes obser- 
vations made on six evenings at wavelength 3880A. ‘The colour of the zodiacal 
light is closely similar to that of the Sun, as has been found by other observers who 
have determined its spectral type. For example, Behr and Siedentopf (3) found 
that it has the colour index of o™-46 at elongation 40°, which may be compared with 
the value of o™-53 for the Sun (10). However, the observations of Behr and 
Siedentopf were made over a comparatively restricted range of wavelengths. 
Elvey and Rudnick (11) deduce from their observations that the colour is very 
close to that of the Sun (although they do not quote an elongation), and a similar 
result has been given by Karimov (12). 

The scan of Fig. 3 shows that the zodiacal light is highly symmetrical about a 
central axis, and our results show that the zodiacal light always has this remarkable 
symmetry when the Sun is not especially active; this result does not agree with 
that found by Behr and Siedentopf, who found marked asymmetries, but it is 
confirmed by the observations of Regener (13). We attribute these discrepancies 
partly to the necessarily greater inclination of the ecliptic to the vertical in the 
work of Behr and Siedentopf, and partly also to disturbance by weak aurorae in the 
higher latitudes. 


TaB_e IV 


Observations of the brightness of the zodiacal light at elongation 40° 
Unit—10~-" times mean brightness of solar disk 


Observer Latitude Height (ft) Date Brightness 
Behr and Siedentopf (3) +46° 11700 1952 3°17 
Barbier (16) +44° 1900 1952/3 3°35 
Roach et al. (7) + 36° 5415 1952/3 5°38 
Regener (9) + 34° 9200 1953/4 5°52 
*Blackwell (17) ee go000 1955 2°90 
Elsasser (14) —29° 4550 1956 2°85 
Divari and Asaad (15) +24° 660 1957 3°26 
Present work (at 6200 A) — 16° 17100 1958 4°42 


(average of all observations) 
Mean brightness at elongation 40° is 3-86 x 1o-“BO. 
*One morning observation only. 


Our measurements are most accurate at elongations in the neighbourhood of 
40°, corresponding to an altitude near 20°, because here the correction for atmos- 
pheric extinction is not great and that for the night sky background is a minimum. 
In Table IV we compare our measurements of brightness at this elongation with 
the results of other investigators, including the results recently obtained by 
Elsisser (14) at the Boyden Station, S. Africa. As the colour of the zodiacal light 
is so close to that of the Sun we have not distinguished between observations made 
at different wavelengths. In reducing other observations we have adopted the 
Stebbins and Kron (10) value of — 26-20 for the photographic magnitude of the 
Sun. There is considerable disagreement between the results of different 
observers even at this favourable elongation, the possible significance of which will 
be discussed in a later paper in this series. 

The corona and the zodiacal light.—The relation between the brightness of the 
zodiacal light and that of the solar corona has been investigated several times. 
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It is now believed that most of the zodiacal light can be regarded as the outer part 
of the solar corona, and this viewpoint is supported by Fig. 7, in which the bright- 
ness of the solar corona and the present measures of the brightness of the zodiacal 
light at wavelength 6200 A are plotted against elongation; the corona measures 
are those of Blackwell (6) made from an open aircraft at a height of 30 000 ft during 
the eclipse of 1954 June 30. The present zodiacal light measures are useful for 
showing this relation, not only because they approach to within 19° of the Sun, but 
also because they were made at the same wavelength as those of the solar corona. 
There remains now a gap of only 5° between the farthest part of the corona that has 
been observed at a total eclipse and the innermost part of the zodiacal light. 
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Fic. 7.—Brightness of the zodiacal light and the corona. 


Measurement of the polarization of the zodiacal light—The measurement of 
polarization is extremely difficult, partly because of the low surface brightness of 
the zodiacal light but chiefly because of the rather uncertain correction for the sky 
background, which increases with increasing elongation. The photoelectric 
measurements of Behr and Siedentopf (3) are outstanding among the early work; 
these measurements were made in the range of elongation 35°<«<go° and 
showed a maximum polarization of about 23 per cent near elongation 35°. The 
range 20° <« < 30° was studied photographically by Blackwell (17) in 1955 using 
an open aircraft flying over the Pacific Ocean at a height of gooo ft and at a position 
where the zodiacal light was vertical at the horizon. Elsiasser (14) has recently 
made new photoelectric measurements at the Boyden Station in the range 
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30°<e<go’, obtaining much larger values of polarization than those found 
by Behr and Siedentopf. 

The measurements at Chacaltaya were made photographically using a grid 
of polaroid strips placed in contact with the photographic plate. The strips 
of the grid were arranged to be approximately radial to the expected position of the 
Sun at the time of a photograph (i.e. about 25° below the horizon), and were cut 
with their axes alternately in a radial and tangential direction so that each strip 
measured the appropriate polarized component of the zodiacal light. A typical 
photograph taken with this arrangement which shows the strong observed polariza- 
tion is reproduced in Plate 3—the actual polarization is two or three times greater 
still than this, the factor depending upon the correction for the sky background. 
Laboratory tests showed that the polarization measures are not unduly affected by 
the large angular aperture of the objective lens. A correction was made for 
deviations of the strips from a radial direction when the Sun is not at a standard 
distance below the horizon. ‘The method is quite sensitive and can easily detect 
polarizations of one or two per cent. Photographs were taken simultaneously 
with the same two filters as were used for the direct photography. 

All of the measurements of polarization for the two wavelengths, corrected for 
night sky background, are shown in an unsmoothed state in Fig. 8 each night being 
represented by at least four or five points, and more when more than one photo- 
graph is available. No significant difference between the polarization values for 
the two wavelengths is apparent. In this diagram we plot also other available 
observations. ‘The present results are in good agreement with those of Blackwell 
(17) and in excellent agreement with the recent results of Elsasser (14), but in 
poor agreement with those of Behr and Siedentopf (3). 
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Fic. 9 —Polarization of the zodiacal light and the solar corona. The corona data are those of 
Blackwell (6) from the 1954 eclipse; the zodiacal light data are an average of the measurements of 
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In Fig. 9 we suggest a working polarization curve formed by combining the 
Chacaltaya values with those of Elsasser and of Blackwell, and we suggest that the 
error of this curve is unlikely to be greater than + 2 per cent of polarization; this 
diagram also shows the relation between the polarization of the zodiacal light and 
that of the corona. 

In Table V we list values of the polarization for an observational model of the 
outer corona and zodiacal light regions based on the data summarized here. This 
model differs greatly from previous models because of the introduction of the new 
data, but we believe that it is more reliable. 


TABLE V 


Polarization of the outer corona (equatorial regions) and the zodiacal light 


Elongation Polarization Elongation Polarization 
4° 3°0 per cent 30° 21°3 per cent 
a 2°8 35° 26°0 
10° 4°0 40° 28:7 
20° 9°7 50° 30'7 
25° 14'8 60° 32°9 
70° 33°2 


The position of the zodiacal light axis—Accurate measurements of the position 
of the axis are very difficultto make. Anaccuracy of atleast + 0°5 is required if the 
measurements are to be useful, and to obtain this accuracy we believe that the 
following conditions of observation are essential. ‘The most important condition 
is that the axis should be nearly vertical at the horizon at the time of observation ; 
if it is inclined to the horizon there is a systematic displacement which depends upon 
atmospheric extinction and is very difficult to estimate. Scarcely less important 
is the condition that the observations are made in a region of low geomagnetic 
latitude where displacements due to aurorae are ataminimum. Also, the observa- 
tions should be made from a high altitude station so that the axis is displaced as little 
as possible from its true position by local variations in extinction, and preferably 
they should be made photographically or photoelectrically. Very few measure- 
ments have been made so far under these conditions. 

The classical observations are those of Hoffmeister (18). ‘These observations 
were made visually from a variety of stations, some of which were within the tropics, 
and the spread of the measurements is about +2°. Thereductions were made very 
carefully, and this probably represents the extreme limit of precision that can be 
obtained visually. The results demonstrated that over a large range of longitude 
the axis deviates significantly from the ecliptic, the maximum deviation being 
about 2°. 

Photographic measurements have been made by Donitch (19, 20) from 
Tamanrasset (in S. Algeria) using an f/4 guided camera and an exposure of 50 
minutes. The observations were made at intervals between 1946 and 1949 and 
showed that the inclination of the axis to the ecliptic is about 2°, and the longitude 
of the ascending node is about 120° (in poor agreement with Hoffmeister who found 
values between 40° and 69°), although both of these quantities seemed variable. 

These data appear to us to be quite convincing, but Fesenkov (21) has recently 
suggested that the axis coincides with the ecliptic. Ina preliminary report on the 
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recent USSR expedition to Aswan, Egypt, he remarks that ‘‘ the deviation of the 
axis of the zodiacal light from the ecliptic must be completely negligible, assuming 
it exists at all’’. This problem is not a trivial or pedantic one. It is very im- 
portant to establish whether or not there really are systematic deviations, for any 
such are good evidence for an extra-terrestrial origin for the zodiacal light. Because 
of this we have used our photographs to obtain measures of the position of the 
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Fic. 10.—Method of measuring microphotometer traces to find the position of the zodiacal light axis. 
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Fic. 11.—Position of the zodiacal light axis in relation to the ecliptic. 
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axis by making microphotometer scans across them parallel to the horizon and 
locating the position of maximum density, allowing for the presence of the brighter 
stars on the scans. The method used for measuring a typical trace is shown in 
Fig. 10, which is self-explanatory. The results of these measurements using only 
the photographs taken through the OB1o (blue) and OY1 (orange) filters, and 
between elongations 30° and 50°, are shown in Fig. 11 for all the dates of observa- 
tion; here the deviation from the ecliptic is plotted against elongation. It will 
be seen that this method locates the position of the axis to within 0°-5 of ecliptic 
latitude and, further, there is no shift of the axis with wavelength. Although the 
range of longitude is not great, there seems to be a significant deviation of the axis 
from the ecliptic. According to our measures, the plane of the zodiacal cloud 
intersects the plane of the ecliptic at longitude 115°, and is inclined at an angle of 
about 1°-5. Both of these quantities are in reasonable agreement with the values 
quoted by Donitch. On Fig. 11 we also plot the position of the invariable plane 
of the solar system, and remark that it is close to the plane of the zodiacal cloud. 
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OBSERVATIONS OF THE ZODIACAL LIGHT FROM A VERY HIGH 
ALTITUDE STATION 


II. ELrecTRON DENSITIES IN INTERPLANETARY SPACE 
D. E. Blackwell and M. F. Ingham 


(Received 1960 October 6) 


Summary 


The evidence for and against the existence of an interplanetary gas is 
reviewed. Spectra of the zodiacal light obtained in 1958 at the high 
altitude station of Chacaltaya (height 17100 ft) in Bolivia with a dispersion 
of 38 A/mm are described and used to investigate the concentration of free 
electrons in the plane of the ecliptic at 1 A.U. from the Sun. No electrons 
could be detected along a line of sight at elongation 32° and it is deduced 
that the electron density at 1 A.U. is not greater than 120 cm=°. 





Introduction.—In Paper I of this series we mentioned the hypothesis that 
the number density of free electrons in the plane of the ecliptic at the distance 
of 1 A.U. from the Sun is about 600cm-*. This hypothesis is now widely 
assumed to be correct, and in the first part of this paper we give a critical review 
of the evidence for and against it. In the second part we describe our observations 
of the zodiacal light spectrum made at Chacaltaya in 1958 and the conclusions 
about interplanetary electron densities that may be drawn from them. 

The evidence for interplanetary electrons.—Some of the evidence has already 
been considered by Blackwell (1), Elsasser (2), Cole (3) and various contributors 
te ~ Geophysical Discussion held in London (4). The first suggestion that the 
eleiron density at 1 A.U. from the Sun may be as great as 600cm~ arose from 
the accretion theory of the solar corona put forward by Hoyle, Bondi and Lyttleton 
(5) in 1947. 

(i) The evidence from the zodiacal light.—In 1953 Behr and Siedentopf (6) 
sought to obtain the electron density far from the Sun using zodiacal light data. 
They then made the first accurate measurements of the degree of polarization 
of the zodiacal light along the ecliptic and found a maximum value of about 23 
per cent near elongation 35°. Until this time it had been universally supposed 
that the zodiacal light is due entirely to the scattering of sunlight by dust contained 
within the solar system, but Behr and Siedentopf suggested that its great polariza- 
tion could not be accounted for by dust scattering alone and that therefore there 
must also be a contribution from electron scattering. ‘The assumption that the 
light scattered by the dust is of zero polarization leads to electron densities of 
the order of 600cm~-* at 1 A.U. from the Sun (6, 7). 

The recent measurements of the polarization of the zodiacal light by Elsasser 
(8) and the present authors (see Paper 1) render the Behr and Siedentopf 
argument even more forceful. These observations give a maximum value of 
33 per cent, which is correspondingly more difficult to explain by dust scattering 
alone than the value of 23 per cent derived by Behr and Siedentopf. Even so, 
the argument for the existence of an electron scattered component still depends 
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entirely upon an assumption that, plausible though it is, has never been tested. 

We stress that acceptance of the Behr and Siedentopf hypothesis inevitably 
implies acceptance of the symmetry properties of the electron cloud. As already 
explained, the zodiacal light is symmetrical about a plane which is very close to 
the ecliptic. If the electron density in this cloud at 1 A.U. from the Sun is as 
great as 600cm~* the electron-scattered component of the zodiacal light contributes 
about one half of the total brightness near elongation 35°, and therefore it also 
must be nearly symmetrical about the ecliptic, implying that the electron cloud is 
symmetrical about the plane of the ecliptic. It may be that such a symmetry 
is unacceptable—and indeed this, as demonstrated later, forms an argument 
against the existence of such a cloud of electrons. The argument can also be 
inverted. If other kinds of observation require the existence of a cloud of 
electrons of this density—as observations of comets apparently do—then we 
deduce that such a cloud must contribute appreciably to the zodiacal light and 
must therefore be nearly symmetrical about the plane of the ecliptic. If the 
original observations do not show that the electron cloud has this symmetry, or 
if the cloud cannot reasonably be expected to show it} then the interpretation in 
terms of a cloud of electrons must be wrong. 

The work of Beckers (g) on the zodiacal light spectrum has apparently con- 
firmed the Behr and Siedentopf hypothesis. Beckers has compared the central 
intensities of the Fraunhofer lines in the zodiacal light spectrum with the inten- 
sities of the same lines in the solar spectrum, using for this purpose a spectrum of 
the zodiacal light obtained by Hoftmeister (10) in 1935. Then, supposing that 
the interplanetary dust scatters the solar spectrum almost unchanged, and that the 
interplanetary electrons scatter a quasi-continuous spectrum in which the 
Fraunhofer lines have been almost obliterated by the large random velocities of 
the electrons, a separation of the dust and electron scattered components may be 
made. According to Beckers the result is that 50 per cent or more of the zodiacal 
light is due to electron scattering. ‘This work will be considered in more detail 
later in this paper, but here we remark only that we believe it to be incorrect 
because proper allowance has not been made for the night sky background to the 
zodiacal light. 

There is evidence, although all of it necessarily indirect, that interplanetary 
dust may indeed scatter highly polarized light. Some of the evidence is provided 
by cometary studies, for some of the light of a comet is due to the scattering of 
solar radiation by dust in its head, and this part may be separated from the 
remainder by a suitable choice of spectral region. Measurements on Comet 
Arend-Roland (1956h) and other comets show a maximum polarization of 
about 25 per cent (11). Admittedly, this is insufficient to account for the polari- 
zation of the zodiacal light, especially as the zodiacal light measurements are 
effectively a weighted average over a large range of scattering angles, but it goes a 
considerable way towards doing so. In any case, good agreement would not be 
expected because the scattering processes and the nature of the scattering material 
may not be the same in a comet head as in the zodiacal light. Fesenkov (12) 
has computed the expected polarization assuming that the interplanetary dust 
is similar to atmospheric aerosols; this gives the value of 28-5 per cent at 
elongation 70°. He infers from this that the Behr and Siedentopf hypothesis 
is an unnecessary one, although it may be correct. A similar conclusion has 
been reached by Whipple and Gossner (13) and by Kloverstrom and Rense 
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(14), on the basis of laboratory measurements, while van de Hulst (15) has 
considered that much, at least, of the polarization may be attributed to dust 
scattering. Opik (16) was one of the first to question the validity of the Behr 
and Siedentopf hypothesis. 

If we except Beckers’ result, which we believe to be wrong, Fesenkov’s 
conclusion fairly assesses the direct observational evidence for an appreciable 
electron component in the zodiacal light. On the other hand, if the electron 
component is assumed to exist we have to explain its symmetry about a plane 
which is so close to the ecliptic. If we suppose that the electrons, with the 
interplanetary gas, are flowing into or out of the Sun it is difficult to see how, in the 
short time available, they can become symmetrical about such a plane. Even if 
we suppose them to be permanent, or semi-permanent, members of the solar 
system, then it is equally difficult to see why the very weak gravitational forces 
associated with the plane of the ecliptic should dominate the random electro- 
magnetic forces which are surely present. These difficulties are most easily 
removed at their source by abandoning the assumption that there exists an 
electron cloud of sufficient density to contribute significantly to the zodiacal 
light. ‘This does not exclude the temporary existence of relatively narrow low 
density corpuscular streams. 

(ii) The evidence from comets.—This evidence is based chiefly upon obser- 
vations of the apparent motions of condensations in comet tails of type I. It 
is commonly supposed that these condensations are moving in orbits under the 
action of a repulsive force which is yz times as great as the gravitational attraction 
of the Sun at this distance, so that analysis of their motion gives the value of 
1—p. ‘There have been many determinations of this quantity, and we mention 
here those of Bobrovnikov (17) for Halley’s comet (1910 II) and the list given 
by Wurm (18); according to Wurm the mean value of 1—y is about 100, 
occasionally it increases to several thousand (17). Biermann (19, 20) has 
suggested that these repulsive forces originate in interaction between corpuscular 
radiation and the plasma forming the comet tail. He finds that the acceleration 
« of a cometary ion exposed to a corpuscular stream is given by 


e2 


“= -neVe, 





omy, 


where mm is the mass of a molecule, me the electron density in the corpuscular 
stream, and a is the electrical conductivity which may be computed from the 
cross-sections for collisions between electrons and molecules, and electrons 
and protons. Inserting the appropriate values, this relation gives the acceleration 
as equal to o-Imecmsec~*. At 1A.U. from the Sun, the acceleration due to 
the gravitational attraction of the Sun is about o-6cmsec~* The mean accelera- 
tion in comet tails is therefore about 60cm sec~*, which corresponds to a electron 
density in the corpuscular stream of about 600cm~*, The agreement with the 
value given by Behr and Siedentopf from the zodiacal light work is remarkable. 

We would criticize this use of cometary data partly on the grounds that it is 
by no means certain that we are observing mass motion in the tail; the apparent 
motions of condensations may only be due to a changing excitation. Even if it is 
demonstrated that these really are mass motions we doubt whether the observed 
repulsive forces may always be attributed to the action of the Sun. The varying 
repulsive forces observed for condensations in different parts of a comet suggest 
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the presence of other mechanisms; for Halley’s comet the condensations of the 
tail showed values of 1 —u between 20 and 150, the average value of 1 — derived 
from the curvature of tail I (the ion tail) was 10 to 20, straight jets near the nucleus 
gave values between 1 and 6, while the envelopes showed approximate balance 
of repulsive and attractive forces. 

(iii) The evidence from magnetic storms.—Biermann (19, 21) has also suggested 
that the observed value of the magnetic activity on quiet days in high geographical 
latitudes may be accounted for by supposing that it is due to a steady corpuscular 
stream of density 10% electrons cm~* at 1A.U. However, we regard this as a 
great over-estimate, for it is based on the Unsdld and Chapman (22) value of 
10°cm~* for the particle density ina great magneticstorm. This value is certainly 
too great, for it, in turn, is based on the work of Briick and Rutllant (23) who 
attempted to observe a displaced K absorption line arising from these particles. 
In his photoelectric work, Smyth (24, 25) was unable to find any increase in 
absorption, which suggests an upper limit of 3 x 10? cm~? for the particle density 
(x). In any case, such a high value as 10°cm~-* is not supported by the 
observations of Paper III. 

(iv) The evidence from the etching of meteorites.—Whipple and Fireman (26) 
have placed an upper limit on the rate at which meteorites have been etched by 
corpuscular radiation in interplanetary space; this limit is 3 x 10~’ cmyear~. 
Supposing that the etching is by protons travelling from the Sun at an average 
speed of 10°cmsec™! (i.e. with energy 5000eV), Whipple (27) finds that the 
corresponding upper limit to the proton density is about 40cm~°. 

(v) The evidence from cosmic radiation.—It has been shown by Brunberg 
and Dattner (28) that the diurnal variation of cosmic ray intensity observed at 
the surface of the Earth is due to an anisotropy of the primary radiation 
reaching the Earth. These authors suggest that this anisotropy is caused by 
the co-rotation of the cosmic ray particles and the solar magnetic field. It is 
now suggested by Block (29) that this, in turn, means that the interplanetary 
gas is also in co-rotation. Block infers from this that the interplanetary gas 
density cannot be greater than o-1cm~* because a gas of density greater than 
this could not be kept in co-rotation by any reasonable solar magnetic field; 
this value of 0-1 cm~* corresponds to an Alfvén field with a strength of 12 gauss 
at the Sun’s poles, which is much too high to be acceptable. It seems to us 
that the weakness of this conclusion lies in the original inference that the cosmic 
radiation is in co-rotation with the Sun. 

(vi) The evidence’ from ionospheric measurements.—In his pioneering work on 
radio whistlers Storey deduced that the electron density at a distance of 3 radii 
from the Earth is about 400cm~* (30). This value however is rather uncertain 
because the method used by Storey gives the total number of electrons along 
the path of the whistler; consequently, the value at the apex of the path depends 
rather critically on the assumed distribution of densities in the lower regions. 
Whilst there have been refinements of the technique by other workers—for 
example Allcock (31) has combined data obtained over a range of geomagnetic 
latitudes—precise knowledge of the electron density in the ionosphere beyond 
about 2 radii depends upon observations of ‘‘nose’’ whistlers, which give 
directly the electron density at the apex of the path (32). Helliwell* finds 
from a study of this type of whistler that the electron density at 5 radii is of 
the order of 1o0ocm~*. 


* R. A. Helliwell, private communication. 
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This short review shows that the electron density in interplanetary space is 
still uncertain by several orders of magnitude, and especially there is a dichotomy 
between the astronomical and the ionospheric data. None of the astronomical 
methods so far described is free of criticism and an independent method of 
investigation is required. We now turn to our own astronomical observations. 

Tne method of measuring electron density used at Chacaltaya.—The method 
is similar to that used several times already for separating the F (dust) and K 
(electron) components of the inner corona. In this application there is inde- 
pendent evidence that the electrons of the corona are at a high kinetic temperature 
so that the light which is scattered by them has a virtually continuous spectrum, 
but the interplanetary dust may be expected to scatter light which, over small 
wavelength regions, has a spectrum identical with that of the Sun. A comparison 
of the central intensities of the Fraunhofer lines in the spectra of the Sun and 
corona therefore gives directly the ratio of the two components of the scattered 
light. The application of this method to the zodiacal light is not so straightfor- 
ward because we do not know the kinetic temperature of the electrons, if they 
exist, at the relevant distance of about o-5 A.U. from the Sun. However, we 
shall show later that the method can be made almost independent of this 
knowledge. 

Assuming now that there are free electrons at 0-5 A.U. from the Sun, we 
suggest that it is unlikely that they have a kinetic temperature of less than 2 x 104 
degrees. This is the value derived by Chamberlain (33) from an application 
of kinetic theory to the interplanetary gas. Values suggested by other authors 
are greater than this. Assuming hydrostatic equilibrium, Pottasch (34) derives 
the temperature of 2-3 x 10° degrees at o-1 A.U., while his curve extrapolates 
to 1:0 x 10° degrees at 0-5 A.U. Chapman (35) has calculated the temperature 
out to 1 A.U., using an empirical function for the variation of angular velocity 
of rotation with distance from the Sun and assuming that the temperature 
gradient beyond 0:26A.U. is close to the adiabatic gradient. His value at 
o5A.U. is 1°2x 10° degrees. 

The temperature of 2 x 10* degrees, the lowest likely value for the kinetic 
temperature, implies that the spectrum of the sunlight scattered by the electrons 
cannot be regarded as completely featureless. As an illustration of the weakening 
effect that may be expected we consider a line which in the solar spectrum is 
observed to have a central intensity of 50 per cent and an apparent width, at 
half intensity, compounded of natural width and instrumental broadening, of 
5A. At the temperature of 2 x 104 degrees and wavelength of 4500 A the whole 
width of the Doppler broadened line at half intensity is 20 A, corresponding to 
an apparent width of about 25A. The depth of the broadened line may be 
computed by equating its area to that of the unbroadened solar line; the result 
is that a broadening of this amount reduces the depth by the factor 9:1, the 
depth of the broadened line being 5-5 per cent. 

We next consider the combination of such a broadened line, which is due 
to electron scattering, with an unbroadened line due to dust scattering. The 
results of computations of the central intensity are given in Table I for a line 
originally of central intensity 50 per cent, for the three kinetic temperatures of 
infinity, 10° and 2x 10* degrees and for three different proportions of the K 
and F components. We deduce from this table that if the ratio K/(K+F) 
is equal to o-5—the approximate value that would be required by the Behr and 
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Siedentopf hypothesis—then the effect of taking the temperature to be 2 x 10! 
degrees instead of infinity is just discernible, while if K/(K+F) is equal to 
o:1 or less the effect is no longer discernible in ordinary photographic photometry. 


TABLE I 


Central intensity of composite Fraunhofer line 


K/(K+F) T=© T=10° T=2 x 104 
o'5 0°25 0°26 0°28 
ol 0°45 0°45 0°46 

° o"50 0°50 0°50 


Alternatively, we may say that if the electron density is great, i.e. 
K/(K+ F)> 0-1, the uncertainty in temperature introduces an uncertainty in the 
density determination of about 10 per cent, but if the density is so low that 
K/(K+F)<o:1 then the errors of photographic photometry are greater than the 
uncertainties of interpretation. If it is our intention only to detect the presence of 
electrons along the line of sight—and it will be shown later that the density is so 
low that only an upper limit can be placed upon it—then the uncertainties in 
interpretation due to ignorance of the temperature are of no account. 

This is an outline only of the method; its difficulties and refinements will 
be explained later. 

The observations.—The spectrum of the zodiacal light has already been 
photographed, chiefly by Fath (36), Slipher (37), Hoffmeister (10) and 
Karimov (38). All of these spectra have been obtained at rather small 
dispersion. 

The spectrograph used in this study has already been described (39). It 
incorporates a Schmidt camera of aperture 8in and focal length 6-3in, the 
correcting plate of which was made by one of us (M.F.I.) at the Cambridge 
Observatories; this, with an 8in x 6in Bausch and Lomb grating of 600 lines 
mm~!, gives a dispersion in the second order at 4300A of 38A/mm. Each 
point in the spectrum is derived from an area of sky of diameter 3-4 degrees. 

Several spectra of the zodiacal light were obtained with projected slit widths 
between 0-13 mm and 0-063 mm, all with the spectrograph raised to an altitude of 
7° above the horizon. The best spectrum for our purpose was photographed on 
103a—O film, specially made by Eastman Kodak, using a projected slit width of 
0-063 mm and an integrated exposure of 7 hours spaced over nights between the 
dates July 4 and July 11; it covered the range 3500-4500 A in the second order. 
As a consequence of the long exposure this spectrum is an average over a range 
of 15° in elongation centred on elongation 30°, the contribution of each element 
in the zodiacal light being weighted according to its intensity. An allowance 
is made for this in the reductions, but as the properties of the zodiacal light vary 
only slowly with elongation this procedure does not introduce a significant 
error. In all of the observations the greatest care was taken to avoid beginning 
the exposure each evening until the last traces of Rayleigh scattering by the 
Earth’s atmosphere had disappeared; in practice, the time of disappearance 
is quite easily determined. The spectrograph was thermally insulated, and 
measurements of a comparison spectrum, put on partly at the beginning and 
partly at the end of the exposure a week later, show that there was inappreciable 
loss of definition through lack of stability. 
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A spectrum of the night sky was obtained subsequently with an exposure of 
7 hours and the same slit width, and a solar spectrum was obtained during 
twilight with an exposure of 30s. All of these photographs were calibrated 
using a standard lamp and a step wedge. ‘The ratio of the zodiacal light bright- 
ness to the night sky brightness was obtained photographically in the manner 
described in Paper I, and an average over the period of observation, weighted 
according to elongation and atmospheric extinction, was used. 

The reduction of the observations.—Although the method is straightforward 
in principle, there are many uncertainties of a practical nature in its application. 
The observed spectrum is a superposition of three spectra of quite different 
origins, namely those of the zodiacal light (appreciably weakened by atmospheric 
absorption), the night sky emission and the twilight emission. The night sky 
spectrum can be observed virtually alone and its contribution to the observed 
spectrum can be allowed for in the usual way by direct photometry. The 
twilight emission spectrum cannot be observed separately and it is best to work 
as far as possible in those regions of the spectrum where its effect seems least 
obtrusive. This means that the D lines of the solar spectrum, which are 
otherwise so suitable, cannot be used in this investigation. ‘The Ha line is 
unsuitable because of the possibility of emission from the night sky and also 
because close to the horizon the night sky spectrum is very complex in this 
region. ‘The H and K lines of Cat+ must also be treated with caution because 
of the occasional presence of a weak emission component of unknown origin (40). 

A photograph of the spectrum of the zodiacal light made with a slit width 
of 0-063 mm is reproduced in Paper I (Plate 3) together with a spectrum of the 
twilight sky. It is clear that the zodiacal light spectrum as observed against the 
night sky contains Fraunhofer lines, but it should not immediately be deduced 
from this, as has been done (41), that the zodiacal light spectrum itself contains 
Fraunhofer lines. The reason for this is that the night sky spectrum also 
shows these lines, particularly H and K of Ca*, the 4227A line of Ca and 
the Balmer series. These lines show strongly on a high resolution spectrogram, 
and a proper allowance must be made for their effect. That it is also important 
to allow for the emission lines of the night sky spectrum is demonstrated by the 
almost complete absence of the line at 4046 A in the composite spectrum, which 
is filled in by the night sky emission line at 4048 A. The presence of an emission 
line at the centre of the K line is clearly shown in the reproduction. 

The spectra have been scanned with the Cambridge Observatories’ micro- 
photometer, and in Fig. 1 we show a microphotometer record of a zodiacal 
light spectrum (superimposed upon the night sky spectrum) made with a 
spectrograph entrance slit width of 1 mm corresponding to a projected slit width 
on the film of 0-063mm. _ In selecting lines for measurement we have avoided 
those which are near to distinct night sky emission lines, and have preferred 
those which are in regions of low night sky emission, where the correction for 
this emission is smallest. 

The records have been reduced by first correcting them for the presence 
of the night sky spectrum. In this we use the results of Paper I to obtain the 
average contribution of the night sky over a broad band of wavelengths, and 
then interpolate to obtain the contribution at each wavelength. The average 
ratio of the brightness of the night sky to the combined brightness of zodiacal 
light and night sky in the region of 4300A was 0°38. 
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The level of the continuum was added to the solar spectra with the help 
of the Utrecht Atlas; this continuum was then transferred to the zodiacal 
light spectra, making clear allowances for the superposed twilight emission 
spectrum. If an undetected emission line is superposed on a Fraunhofer line, 
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Fic. 1.—Microphotometer record of the zodiacal light spectrum (superimposed on the night 
spectrum). Width of spectrograph entrance slit=1 mm. 








th 


asi 
abs 
eva 


wh 
eqt 


is 1 


isa 
deg 
of | 
thre 
by 1 
of 1 





[22 


elp 
cal 
ion 
ne, 





No. 2, 1961 The zodiacal light from a very high altitude: II 137 


the effect will be an apparent weakening of the Fraunhofer line which will 
simulate the presence of electron scattering. The line depths were then 
measured directly on each of the spectra with the results given in the third 
column of Table II. The results show approximate equality of the line depths 
in the zodiacal light spectrum and the solar spectrum. 














TaBLe II 
Comparison of intensities of Fraunhofer lines in the zodiacal light spectrum and the solar 
spectrum 
A+ AA 
| I,dX Line depth in z.l. spectrum 
: Ao—AA Zodiacal light 
Line ayy - - 
| I,dn Line depth in solar spectrum 
Ao— AA Sun 
Ca 4227 0°96 
Fe 4384 0°94 0°94 
Hy 4341 0°93 1°10 
Fe 4326 IIo 1°20 
*CatH 3968 0°97 
G-band 4308 1°05 I‘Io 
H8 4102 085 0°87 
Average 0°97 1°04 


* Corrected for presence of emission component. 


A disadvantage of this method of reduction is that the result depends upon 
a single measurement of the microphotometer trace at the position of maximum 
absorption. To reduce the uncertainty of the measurement we have also 


evaluated the quantity 
AotAA 
| Ian 


Ao— AA 


rot AA 
| Ian 
ss" 


where A, is the wavelength of maximum absorption and Ad is about 3A. The 
equivalent width of a line, defined as 


| * 1d 
0 


is unaltered by electron scattering, but if AA is not too great a fraction of the 
Doppler width of the broadened line the quantity 

Ag+ La 

| 1,aa 

Ay— AA 
isa good measure of the central intensity. For an electron temperature of 2 x 104 
degrees, AA may be as much as 3 A without introducing an appreciable uncertainty 
of interpretation, while the accuracy of measurement is considerably increased 
through the use of a greater length of microphotometer record. Results obtained 
by this method are given in column 2 of Table II. Combining the two methods 


of measurement, the mean result is 
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The uncertainty in this result is difficult to estimate; we simply suggest that it 
is very unlikely to be in error by more than 10 per cent. This estimate is a 
combination of the errors of photographic photometry with the uncertainty 
in the separation of observed zodiacal light intensity into a true zodiacal light 
intensity and a night sky component. 

Interpretation in terms of the electron density in interplanetary space.—It is clear 
that these observations have not shown the presence of any free electrons along the 
line of sight. We now compute the electron density at 1 A.U. corresponding to 
K/(K+F)=o'1, this being the least electron density detectable by these 
observations. 


Electron [@) 


€ 
Sun Observer 





Fic. 2.— Scattering of sunlight by interplanetary electrons. 


In the notation of Fig. 2, and defining surface brightness as the rate of 
flow uf energy through unit area at the Earth’s surface and contained within 
unit solid angle, we calculate the ratio, ®, of the brightness of the electron cloud 
at elongation ¢« to that of the Sun to be given by 

( ‘7 
e+ | ne(r){x +0086} dé 


167sine ], 
neglecting the finite size of the Sun. Here, o is the Thomson scattering 
coefficient, Q is the solid angle subtended by the solar disk and me(r) is the 
electron density at distance r from the Sun. The value of the electron density 
at a particular distance corresponding to a given value of ® does not depend 
greatly upon the form of me(r); as our final result is that me is small at distances 
of the order of 1 A.U., me varies at least as rapidly as 1/r? and we take it to be 
proportional to 1/r?._ In the computations we use the value for the brightness 
of the zodiacal light of Table IV (Paper 1), i.e. the brightness at elongation 
40° is 3°86 x 10 B., and correct it to the appropriate wavelength using the 
data of Fig. 6 (Paper 1). Putting K/(K+/)<o-1, we find the final result 
me (at 1A.U.)<116cem~, 

This upper limit for the electron density depends upon the assumption 
that the kinetic temperature is not less than 2 x 104 degrees. It seems to us 
that such an assumption is reasonable, but electrons in the concentrations 
required by the Behr and Siedentopf hypothesis could still be detected at 
much lower temperatures through the added presence of wings or line 
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broadening. In Fig. 3 we compare the observed profile of the Hy line at 
4341 A in the solar spectrum with that of the same line in the zodiacal light 
spectrum (corrected for background), and it is evident that there are no 
perceptible wings or a broadening of any kind. The difference in total width 
at half intensity between the two is less than 0-28 A, corresponding to a kinetic 
temperature of only 4 degrees. As this temperature is unreasonably low we 
conclude that in no circumstances are the present results compatible with the 
Behr and Siedentopf hypothesis, and that it is very unlikely that the electron 
density exceeds 120cm~*, 
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Fic. 3.—Comparison of the observed profile of the Hy line at 4341A in the zodiacal light 

spectrum with that of the same line in the solar spectrum. 

The narrow slit observations were made between July 4 and July 11 when 
the average enhancement in brightness of the zodiacal light over this period was 
about 8 per cent, but this change in brightness has not been detected as a change 
in Fraunhofer line intensity. Before changes of this kind can be detected 
spectroscopically it will be necessary to develop faster and more accurate 
methods of observation. 

As already mentioned, our result is a direct contradiction of that obtained 
by Beckers (g). We attribute this chiefly to Beckers’ neglect of the night sky 
background. ‘The ratio of zodiacal light brightness to night sky background 
depends upon the quality of the sky and the conditions of observation, and 
none of these data seems to have been available to Beckers. Also, we regard 
it as essential to make observations of the solar spectrum and the night sky 
spectrum with the same spectrograph as was used for the zodiacal light 
observations, using also the same slit width. 
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There are two outstanding difficulties connected with our result. The first 
is that we must now attribute the whole of the polarization of the zodiacal light 
to scattering by dust particles in interplanetary space. Yet, as shown in Paper 
IV, this polarization is much greater than can be explained by laboratory experi- 
ments with likely scattering materials. 

The second difficulty is related to the work of Biermann and Schliiter 
on the accelerations of condensations in comet tails, which has already been 
discussed. ‘These two authors attribute the accelerations to the effect of a 
steady corpuscular stream of electron density about 600cm~* at 1 A.U. from the 
Sun. While we do not deny the occasional existence of such a corpuscular 
stream we suggest that the accelerations usually observed have a different origin. 

An important consequence of our result is that the reliability of existing 
measurements of electron density in the outer corona is now very questionable. 
These measurements (42, 43) are based upon observations of polarization 
and assume that the polarization of the F (dust) component is negligible in 
comparison with the total polarization. But the total polarization at elongation 
5° is only 2-7 per cent, and if the polarization of the light scattered by dust at 
elongation 35° is 33 per cent, the polarization of the F component at 5° elongation 
may well be as much as 2-7 per cent. Our conclusion is that we have no reliable 
quantitative knowledge of electron densities beyond about 8 solar radii, and 
our only evidence for the existence of electrons beyond this distance is provided 
by the work of Hewish (44, 45) and Vitkevitch (46) on the scattering of radio 
emission from the Crab nebula by the outer corona. However, such scattering 
observations only show the presence of electrons and do not enable a measure- 
ment of density to be made. 
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OBSERVATIONS OF THE ZODIACAL LIGHT FROM A VERY HIGH 
ALTITUDE STATION 


III. Tue DisturBeED ZOpDIACAL LIGHT AND CORPUSCULAR RADIATION 
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Summary 


Photometric observations of the zodiacal light made at Chacaltaya in 
the Bolivian Andes during 1958 show that there are changes in brightness 
which are correlated with geomagnetic activity. The changes in the zodiacal 
light and the surrounding sky during the intense magnetic storm of 1958 
July 8-9 are described. The increase in sky brightness that was then 
observed is tentatively ascribed to scattering by free electrons in a corpuscular 
stream, the electron density in the stream being about 300cm~*. There 
was also a corresponding increase in zodiacal light brightness; the hypothesis 
that this may be attributed to fluorescence of the interplanetary dust on 
impact with the protons of this corpuscular stream is examined, but found 
to be inadequate. 





Introduction.—Whilst at Chacaltaya the 1958 Cambridge Observatories 
expedition* made a special study of the variation with time of the brightness of 
the zodiacal light and the surrounding night sky. Previous studies of this problem 
by other observers had given conflicting data and we hoped that the exceptional 
atmospheric conditions at Chacaltaya would enable us to decide whether or not 
the zodiacal light itself varies in brightness or position. Whilst this was not the 
primary problem, we also wished to examine the possibility of detecting corpuscular 
radiation directly, during its passage between the Sun and Earth, through the 
scattering of sunlight by free electrons in the corpuscular stream. The time of 
the expedition was particularly favourable for studies of this kind, for the observa- 
tions were made at the peak of the solar cycle during a period when the solar 
activity was exceptionally great. 

Previous observations of changes in the brightness of the zodiacal light.—There 
have been many reports of variations in the brightness of the zodiacal light. We 
list some of these below, but at the same time we emphasize most strongly that 
studies of possible variations in the zodiacal light can be made only with great 
difficulty and any results should be treated with extreme caution. Worthwhile 
data can scarcely be obtained except from good photographic or photoelectric 
photometry carried out under excellent meteorological conditions in the tropics, 
and for this reason we reject almost all of the available visual observations, even 


_ * The circumstances of the Cambridge Observatories expedition to Chacaltaya and the observa- 
tions of the brightness, polarization and position of the zodiacal light made there are described in 
Paper I of this series. 
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though they have been made with the greatest care. The principal changes that 
have been reported are as follows. 


(i) Short period fluctations.—Reports of changes occurring during one evening 
within a time interval of afew minutes arecommon. The best known study of this 
kind is that of Jones (1) made between 1853 April 2, and 1855 April 21, during a 
voyage on the Pacific Ocean chiefly between latitudes + 20° and +40°. Jones’s 
observations were made visually whilst he was on board ship, and he recorded the 
brightness and position of the zodiacal light together with changes in its brightness 
which he called ‘‘ pulsations’’. Hulburt (2) has drawn attention to the correla- 
tion that exists between these times of zodiacal light variability and the times of 
magnetic activity as observed at Greenwich. Jones himself was unaware of any 
such relationship, and it is difficult to escape the conclusion that the phenomena 
described by him did actually occur*. There are other reports in the literature— 
see, for example, Hulburt (2)—of changes in the zodiacal light during periods of 
auroral activity, but we postpone until later the question of whether the observer 
is deceived by the presence of aurorae into thinking that the zodiacal light has 
changed. 


(ii) Occasional enhancements of the zodiacal light during one whole night.— 
Visual observations of this kind are not uncommon (2), and we refer here only 
to the unusually bright zodiacal light of 1896 March 4 (3), which occurred during 
a period of auroral activity. Apparent changes of brightness from night to night 
occur in most detailed lists of photometric observations—see, for example, Roach 
(4) and Huruhata (5)—although an opinion is rarely expressed about the reality 
of these apparent changes. 


(iii) An annual variation.—An annual variation of brightness over a range of 
two to one has been observed photoelectrically by Elvey and Roach (6). ‘Thom 
(7) has assembled many visual observations from a variety of sources and has 
found a similar variation. 


(iv) Long period variations.—Huruhata (5) has made photoelectric observa- 
tions over a period of four years, finding an irregular variation from year to year 
amounting at most toa factor ofseven. Thom (7), using the previously mentioned 
visual data, also finds a variation and some correlation with solar activity. Table 
IV of Paper I, which lists absolute brightness values obtained by various observers 
during the period 1952-8, using non-visual techniques, apparently supports this 
idea that there is a variation from year to year. 

All of this evidence is contradicted by Regener (8) who cannot find any signifi- 
cant variation during fourteen months continuous observation over the period 
1953-4. Itis perhaps significant, however, that his observations were made during 
a period of very low solar activity. A daily record of the appearance of the zodiacal 
light at the Harvard College Observatory (latitude 42°N) ever the period 1876- 
1890, has been analysed by Searle (9) who states that the variations which seem 
to have occurred cannot be established with certainty, and any support for a 
correlation between the variation and the occurrence of aurorae is certainly very 
feeble. Barbier (10) also, working photoelectrically at Haute- Provence (latitude 


*He noted in his diary for 1854 January 30, “There can be no doubt that there are pulsations in 
the zodiacal light. 1 noticed them last evening; but, it being Sunday, made no particular record 
of them.”’ 
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+ 44°), finds no evidence for a variation of zodiacal light brightness. He worked 
during the periods 1951-2 and 1952-3 when solar activity was low. 

We conclude that it is not possible to deduce from the existing data whether 
or not the zodiacal light is of constant brightness. 

The Chacaltaya observations.—The requirements for a site suitable for an 
investigation of possible variations of the zodiacal light are very stringent. 
Evidently the site must be near the geographical equator so that the varying 
inclination of the ecliptic to the horizon does not lead to systematic errors, and it 
must also be at a high altitude so that observations may be made close to the 
horizon, and have an excellent climate in which the extinction is low and varies 
little from night to night. In addition, it must be near the geomagnetic equator 
so that there is as little disturbance as possible from aurorae. Chacaltaya satisfies 
all of these requirements. In particular, it is at a low geomagnetic latitude 
(—3°) and so is especially suitable for the attempted observation of corpuscular 
radiation mentioned at the beginning of this paper. 

The method of observing the zodiacal light and the results have already been 
described in Paper I. In this section we attempt to ascertain whether or not the 
observed variation in these data from day to day can be regarded as real. 

Correlation between brightness measurements at \4470A and \6200A.—The 
measurements of brightness in these two spectral regions have been made almost 
completely independently, and we may therefore test the reliability of these 
measurements by examining the correlation between them. To do this we have, 
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Fic. 1.—Correlation between mean surface brightness of the zodiacal light at the two wavelengths 
4500A and 6200A. 
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for each photograph, averaged over all elongations the deviation from the mean 
brightness. In Fig. 1 we plot one set of measurements against the other and 
obtain a reasonably good correlation with a slope of unity, showing that the same 
relative change of brightness occurs in the two spectral regions. 

Further, supposing that a perfect correlation may be expected to exist between 
the two quantities, the diagram shows that the probable error of such brightness 
measurements, averaged over one plate, is about 8 per cent. 

Correlation between zodiacal light brightness and solar activity.—As it now seems 
likely that the observed variations of brightness are real, we attempt to correlate 
them with solar activity, as indicated by the magnetic index K,. We have been 
particularly fortunate during the expedition to have observed the zodiacal light 
before, during and after the magnetic storm following an important solar flare. 
The flare occurred on July 7 at 00"39™ U.T. and was recorded at Mitaka and 
Hawaii Observatories (11) as of importance 3+. The intense magnetic storm 
which followed it on July 8—g was one of the most important during the three year 
period 1957—9, the value of the planetary magnetic index K,, (11) remaining 
at 9 for 15 hours. 

In Figs. 2 and 3 we reproduce the plot of surface brightness against elongation 
for the two colours, which has already been given in Paper I, and show the points 
appropriate to days when the magnetic index K,, was greater than 4 (i.e. exception- 
ally high) at the time of observation; evidently there is a fair correlation between 
magnetic index and brightness in the two colours. 
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Fic. 2.—Surface brightness of zodiacal light at wavelength 6200A showing values appropriate to 
active days. 
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Fic. 3.—Surface brightness of zodiacal light at wavelength 4500A showing values appropriate 
to active days. 
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Fic. 4.—Change of average brightness of zodiacal light as a function of magnetic index. 
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To investigate this relation further we have formed for each day the deviation 
of brightness from the mean averaged over elongation and over the two colours. 
In Fig. 4 these deviations are plotted as a function of magnetic index for the time of 
observation, showing a remarkably good correlation, the total change in brightness 
when the magnetic index changes between o and g being about 40 percent. The 
correlation is not substantially altered when a phase difference of up to one day 
is introduced. These results apparently justify the suspicion that has existed for 
more than a century, that the zodiacal light itself does change during a time of 
auroral activity. We emphasize that this observed change in zodiacal light 
brightness is quite distinct from a change in general sky brightness which may or 
may not be due to a localized aurora*. However, no rapid changes of the kind 
described by Jones were observed, either visually or photographically, and we are 
inclined to doubt whether these changes are associated with the zodiacal light 
itself, remembering that Jones often observed at rather high geomagnetic 
latitudes. 

We have examined other published photometric data for the effect, but have 
found that the data are too scanty and the range of magnetic index too small to 
show a good correlation. 
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Fic. 5.—Relation between brightness of zodiacal light at «=40° and date of observation; 
a correction has been made for magnetic activity using the data of Fig. 4. 


Examination of the data for an annual variation of brightness.—We are now 
able to use the graph of Fig. 4 to allow for the effect of magnetic activity on the 
zodiacal light brightness and to reduce the data to zero magnetic activity; the data 
may then be examined for an annual variation. Our observations are restricted 
to the period June 4—August 2, but this interval should be sufficient to show 
part of an annual variation, if it exists with an amplitude equal to that already 
reported. In Fig. 5 we plot against date the relative brightness of the zodiacal 
light at elongation 40° reduced to zero magnetic activity—this elongation being 
chosen because the brightness values are most accurate here. Each point on the 


* The difficulty of distinguishing between a change in zodiacal light brightness and a change in 
sky brightness is further discussed on p. 150. 
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graph is a weighted mean of observations made in the blue and red regions. All 
of the photographs were taken at approximately the same time in the evening and 
reduced in the same way; a correction was made for extinction, using for this 
purpose observations of stars made with a photoelectric telescope. The graph 
shows no evidence for an annual variation of brightness such as has been observed 
by Elvey and Roach (6). 

The position of the symmetry axis of the disturbed zodiacal light.—Six repre- 
sentative points on the sy’ metry axis for the dates July 8, 9 (when the magnetic 
index was greatest and the zodiacal light most disturbed) are shown plotted in 
Fig. 6, together with data from all other dates; the points for the two active days 
are averages for red and blue light. The scatter of all the points indicate the 
probable error of the data obtained on inactive days, and we conclude that the 
zodiacal light was not significantly displaced from it normal position on these two 
occasions. 
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Fic. 6.—Symmetry axis of quiet zodiacal light and points obtained on active days. 


Although we do not give isophotes for the zodiacal light, our scans show that 
there was not a significant change in the shape of these isophotes on the active 
days, i.e. the zodiacal light cone brightened as a whole. 

The brightness of the night sky.—We quote here the night sky brightness at an 
altitude of 15° above the horizon, both on the ecliptic and at about 45° on each 
side of the ecliptic. 'The greatest change of brightness during the magnetic storm 
of July 8—g was observed north of the ecliptic. To illustrate this we plot in Fig. 7 
as a function of date the relative sky brightness measured in the red and blue 
regions, together with the planetary magnetic index K,. There is evidently a 
close correlation between the brightness data for the two wavelength regions, 
and also a correlation between brightness and magnetic index; the increase 
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in brightness at the time of the storm on July 8~9 is particularly striking. The 
background sky brightness on the ecliptic, obtained by interpolation between the 
background at the sides of the zodiacal light, also changed, but it shows a much 
weaker correlation with magnetic index, the sharp rise on July 7 being scarcely 
discernible. The same is true of the brightness of the sky south of the ecliptic. 
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Fic. 7.—Relation between sky brightness at 4500A and 6200A north of the ecliptic, 
and the magnetic index. 


The following diagram, Fig. 8, summarizes these observations of the effect of a 
solar flare on the zodiacal light and the neighbouring night sky. The diagram 
shows typical scans across the zodiacal light made under quiet conditions and 
soon after a solar flare; the level of the night sky background is also shown. 
The brightness of the zodiacal light itself increased but its position remained 
unchanged ; these effects are superimposed on a changed night sky background. 
It is clear from this diagram that it is difficult to obtain the true brightness of the 
zodiacal light under disturbed conditions because of the uncertainty of the inter- 
polation of the sky background between measures made north and south of the 
ecliptic. In our work we have made a linear interpolation. Because of this 
difficulty all observations of this nature should be treated with great caution. 
The most reliable data are those which refer to small elongations where the ratio 
of zodiacal light to the sky background is a maximum. 

The colour of the increment of light from the northern night sky.—This is difficult 
to establish with precision. Using the data of Fig. 7 and allowing for extinction, 
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Fic. 8.— Effect of a solar flare on the zodiacal light and night sky brightness. The diagram shows 
representative scans across the zodiacal light before, and soon after, a flare. 
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Fic. 9.—Position of solar flare on July 7 in relation to the evening sky. 














152 D. E. Blackwell and M. F. Ingham Vol. 122 


we find that the extra light which comes from the northern night sky during 
enhanced magnetic activity has approximately the same colour as the solar disk. 
The actual ratio observed is 
Brightness of increment of night sky at o-62u /Brightness of Sun at 0-62 
Brightness of increment of night sky at 0-45u/ Brightness of Sun at 0-454 





= 1°15 + 0°20. 


The large probable error is a measure of the difficulty of this observation. 

Interpretation of the observed effects of the solar flare.—We begin by describing 
the circumstances of the flare of July 7. The flare occurred at the following 
position on the solar disk —N24°, Wog°. Examination of Fig. 9 shows that a 
jet of corpuscular radiation emitted from the region of this flare in a direction normal 
to the solar surface would be directed away from the Earth, and carried away 
further by the solar rotation. Nevertheless, the occurrence of an intense geo- 
magnetic storm shows that the Earth must have passed through at least the outer 
parts of the stream during July 8-9. The corpuscular stream was also responsible 
for marked perturbations in the motion of the Satellite 195851 (12). The 
increased brightness of the northern sky observed in the period June 9 onwards 
(see Fig. 7) is probably also due to the same active group, although no important 
flares were recorded during this period. 

The observed increase in general sky brightness may be attributed either to an 
aurora or to the scattering of sunlight from the free electrons contained in a 
corpuscular stream. We incline to the second explanation for the following 
three reasons. First, the increase in sky brightness is observed chiefly in the part 
of the sky which is north of the ecliptic and this is compatible with scattering from 
a corpuscular stream emitted from an active region on the northern hemisphere of 
the Sun. Second, the colour of the additional sky light is close to that of the Sun, 
whereas aurorae often appear red or green. ‘Third, aurorae occur only very 
rarely at the low geomagnetic latitude of 3°S. Alternatively, it might be argued 
that if the extra sky brightness is really due to scattering in interplanetary space, 
the increase should have started soon after the flare itself—whereas the brightness 
increased with the magnetic activity. This argument is not correct. According 
to Fig. 9, the corpuscular stream was initially directed away from the Earth, and the 
diagram of Fig. 10, which shows the corpuscular stream projected on to the plane 
of the ecliptic, shows that the effect need not be expected until at least one day 
after the flare. ; 

Supposing that the effect is really due to scattering by corpuscular emission, 
we may calculate from the increase in brightness the order of magnitude of the 
electron density in such a temporary corpuscular stream, which should be distin- 
guished from a more or less permanent interplanetary gas. If, on the other hand, 
the effect is not wholly due to electron scattering, the density obtained represents 
an upper limit. Taking the data already given, and assuming that we are looking 
through a total path length of about o-5 A.U. in the corpuscular stream, the density 
would be about 300 electrons cm~*, assuming that it is independent of distance from 
the Sun. Such a value is not incompatible with those already suggested on other 
grounds. For example, Biermann (13) has suggested that the disruption of 
comet tails is caused by solar corpuscular radiation, and he supposes a density 
between 10% and 10°atomscm~*, Chamberlain (14), using the observed intensity 
of Ha in aurorae, has proposed the value of 1atomcm~*. For a magnetic storm 
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showing a maximum value of the magnetic index K, equal to g, Ferraro (15) 
deduces the atom density in a cloud of corpuscules to be 40cm~%, although 
elsewhere he suggests (16) an upper limit of 1oocm~*. 

We suggest that these two effects which can contribute to an increase in the 
sky brightness after a flare, i.e. an aurora or scattering from temporary corpuscular 
radiation, may be distinguished by the change in the polarization of the sky light. 
If the extra light is auroral in origin, it would be unpolarized ; whereas if it is due to 
electron scattering, it would be strongly polarized. Of course, such a separation 
is only practicable in tropical regions*, where aurorae are very infrequent. 


Corpuscular 
stream 








— Earth Sun 


Fic. 10.—Geometry of the emission of a corpuscular stream on Fuly 7. Projection on to the plane 
of the ecliptic. 


The increase in the brightness of the zodiacal light.—It seems significant that 
when the zodiacal light increases in brightness its position in space remains 
unchanged. ‘The position of the zodiacal light cone, between the ecliptic and the 
invariable plane of the solar system, suggests that its source is a permanent feature 
of the solar system and is primarily gravitationally controlled, i.e. that free electrons 
are not an important contribution. This is supported by our spectroscopic 
measurements. If the light cone does not move when it brightens, following a 
solar disturbance, it seems unlikely that the extra light comes by electron scattering. 
We now examine the suggestion that this extra light is due to a fluorescence of the 
dust particles caused by collision between the high speed protons of the corpuscular 
stream and the dust particles of the zodiacal cloud. 

To investigate the possibility that such a mechanism is of importance, we 
calculate the expected increase in brightness of the zodiacal light assuming that 
when a proton contained in the corpuscular stream collides with a dust particle 
a certain fraction of its kinetic energy is converted into energy of radiation in the 
visible spectrum. We neglect the electrons because of their small kinetic energy. 
Because of the various uncertainties, the calculations are necessarily very crude. 


* This method of separation will be attempted during an Oxford—Cambridge expedition to 
Chacaltaya in 1961. 
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In the calculation we assume that the corpuscular emission has a density of 
300 protons cm~*, which is independent of distance from the Sun, corresponding 
to our estimate of an electron density of 300cm-%, and that the corpuscles are 
travelling at a speed of 4x 10°cmsec™1. Such a velocity as this has been sug- 
gested by Meinel (17) on the basis of studies of the broadened H« line in the spectra 
of aurorae. 

We also need the total cross section A, of all dust particles in unit volume of 
interplanetary space. In the calculation of this we anticipate a result from Paper 
IV of this series, where it is supposed that the distribution of particle size is given 
by the relation, 


n(a)da=Ca~? da 


where n(a)da is the number of dust particles having radius between a and a + da. 
We suppose also, that the space density of dust particles varies with distance from 
the Sun as 7r-*. 

Following the methods of Paper IV, we find that proton impact on the dust 
cloud gives an increment in brightness at elongation « equal to 


6, Xe 
1°3 x 10° 7Rf cosec**1 | sin*—* 6 dé | Ca*-P da_ erg sec~! cm~ sterad- 
6, a 


where f is the fraction of the kinetic energy of the protons which is converted into 
radiation in the visible spectrum on collision with a dust particle. Evaluation of 
this expression, with C=4:28x10-%, R=1-5x10!%cm, a=4x 10cm, 
a=1, p=5, leads to an increment which is only f/100 of the brightness of the 
zodiacal light at elongation 35°. If {<1 we conclude that this mechanism is not 
able to explain the observed increment. 

We are faced here with an apparent contradiction. After a solar flare there is 
apparently a brightening of the zodiacal light which must surely be attributed to 
increased emission from the zodiacal dust cloud. If the increased emission is due 
to fluorescence, we require a proton density which is greater than the maximum 
electron density permitted by the changes in the night sky brightness. A way out 
of this difficulty could be to suppose that the corpuscular stream contains a high 
proportion of neutral hydrogen which does not scatter solar radiation, but which 
can excite fluorescence. The calculations of Kahn (18) of the ionization of 
hydrogen in a corpuscular stream do not support this view. Yet, as we have 
mentioned already, Biermann has suggested that the observed accelerations of 
ions in comet tails can be accounted for only by assuming a density of corpuscular 
radiation of between 10° and 10° cm-*, the latter corresponding to an intense 
storm. We do not believe that such high densities of free electrons can exist in 
interplanetary space because unless they occur in very thin clouds the resulting 
increase in sky brightness could hardly escape detection. If Biermann’s hypo- 
thesis is correct a corpuscular stream must contain a high proportion of neutral 
hydrogen. This is just the condition that we require to explain the increase in 
zodiacal light brightness through fluorescence without a large increase in sky 
background through electron scattering. 

Conclusions.—We find that following a solar flare there is an increase in zodiacal 
light brightness and also an increase in the brightness of part of the sky. The 
increase in sky brightness is attributed to the temporary emission from the Sun 
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of a beam of corpuscular radiation of density 300 electrons cm~*, and this view is 
supported by the colour of the increment in sky brightness and the place in the sky 
where the increase occurs. ‘The increase can scarcely be due to an aurora because 
of its colour and because of the rarity of aurorae at the low geomagnetic latitude 
of 3°S. 

The increase in brightness of the zodiacal light is not accompanied by a shift 
in its position, which suggests that the extra emission is from the interplanetary 
dust cloud. This extra emission cannot be due to interaction between the dust 
cloud and the protons of the corpuscular stream because these carry insufficient 
energy. While it is possible that the increase is due to fluorescence of the inter- 
planetary dust under the action of short wavelength radiation, the scarcity of the 
data precludes adequate discussion of this hypothesis. 


The Observatories, 
Madingley Road, 
Cambridge : 
1960 October 5. 
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OBSERVATIONS OF THE ZODIACAL LIGHT FROM A VERY HIGH 
ALTITUDE STATION 


IV. THe NATURE AND DISTRIBUTION OF THE INTERPLANETARY DusT 
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Summary 


The problem of the separation of the K-(electron) and /F-(dust) 
components of the corona is reconsidered in the light of the revised value 
of the interplanetary electron density at 1 A.U. derived in Paper II of this 
series. A method of separation is developed in which it is not necessary 
to assume that the polarization of the F-corona is zero, but instead a model 
of the electron density is assumed. A model of the surface brightness and 
polarization of the corona and the zodiacal light is constructed on the basis 
of existing observations and the method of separation applied to it. The 
results are compared with theoretical expressions for the surface brightness 
of the light scattered by solid particles by (a) diffraction and (5) reflection 
and refraction. Values of the albedo, mass density and number density of 
the particles composing the zodiacal cloud are derived. The theory does not, 
however, account for the observed reddening of the corona. 





1. Introduction.—In Paper II of this series (1) it was shown that the density of 
free electrons in interplanetary space at 1 A.U. from the Sun is at most 120cm~* 
and probably considerably less than this. It is therefore certain that the inter- 
planetary gas cannot account for the whole of the observed polarization of the 
zodiacal light and probably for not even a significant fraction of it. Hence the 
greater part, if not all, of the polarization must be due to scattering by interplane- 
tary dust. Hitherto, computations of the electron density in the outer corona 
(r>2R,,) have been made on the assumption that the polarization of the F-corona 
isnegligibly small. Now, however, this assumption has been called into question 
since, although at elongations less than 5 R,, the F-corona is probably due entirely to 
diffraction and has zero polarization, at larger elongations the strongly polarized 
component due to reflection and refraction becomes increasingly important. 
Furthermore, if the estimate of the electron density at 1 A.U. is to be greatly 
reduced, a revision of this quantity in the coronal regions near the Sun becomes 
necessary. ‘The purpose of this paper is to develop and apply a method whereby 
the corona and the zodiacal light may be separated into their two components in the 
light of the revised values of the electron density at 1 A.U. and without making 
any assumptions about the polarization of the F-corona. 

2. Outline of the method.—lf we denote by K and F respectively the surface 
brightness of the K- and F-components of the corona, by p, and p, their polariza- 
tions and by pf, the observed polarization of the total corona, then we have the 
simple relation 

Po= poet Fer (1) 





in which the quantities which are observed arep,andK+F. Ifweputp,=othen 
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it is possible to deduce both K and p, by a method of successive approximation 
due to van de Hulst (2). If, however, p,#0 this method is no longer available 
and we have therefore used a different means of separating the two components 
based on the following considerations. 

Close to the Sun (at elongations < 2R_) it is possible to perform the separation 
by spectroscopic means without appeal to assumptions concerning polarization. 
The same situation regarding the separation pertains, as was shown in Paper II 
of this series, at great distances from the Sun in the region of the zodiacal light, 
whose spectrum has now provided an upper bound to the electron density. Thus 
close to the Sun (i.e. for 1<r<4R,) it is in principle possible to obtain fairly 
reliable values of n,(r) and at great distances there is a useful upper bound. 

We now construct a series of functions n,(r) which join on smoothly to the 
values in the range 1 <r<4R_ and whose values at r=1 A.U. are less than the 
upper bound. A given function n,(r) leads to values for K and p,. Knowing 
(K+F), F is obtained by subtraction and, from equation (1), p,/. Division 
yields p,, and the separation is achieved. The mathematical details of this method 
and the criteria for deciding which of the functions 7,(r) leads to the correct forms 
of K, F, p, and p, will be considered in the following sections. 








Fic. 1.—Geometry of scattering in the corona and interplanetary space. 


3. Basic equations for the scattering of sunlight by small particles in interplanetary 
space.—The analysis presented here is based on the work of van de Hulst (3). 
We define the scattering function 1(@) of a spherical particle as the amount of 
radiation scattered per unit solid angle in a direction inclined at an angle @ to the 
direction of incidence, divided by the incident radiation that hits the geometrical 
cross section of the particle. Let.J, dA denote the flux of the solar radiation in the 
wavelength range A to A+ dA at a distance R=1A.U. from the Sun. Consider a 
column, lying in the line of sight, of height dA, cross section a, seen at elongation «. 
Let its distance from the earth be A and from the Sun be r (Fig. 1). 

Further, let there be n(a)da particles per cubic centimetre at 1 A.U. whose 
radius lies between aanda+da. These particles receive solar energy at a rate of 


2 a 
J, a(*) oddn(a)da( na* erg sec 


where we have assumed that the space density of the particles varies with distance 
from the Sunasr-*. Of this energy a fraction 1(@)A~ is scattered to 1 cm? of the 
surface of the earth. This energy is seen coming from a solid angle cA-* so that 
the surface brightness of the column is 


2+a 
Fi4=d, a(*) n(a)dana*I(@)dA erg cm~* sec~ steradian-, 
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Now from Fig. 1 


sine sine 
nia and nhs 


and so for the surface brightness at elongation « we have 


dé 








F,(c)=J ,dAR cosec!+%e | ‘a | "n(a)ra® sin* 61(0) dO da (2) 


where @,, a, are the radii of the smallest and largest particles present. The 
expression on the right of equation (2) may be divided into two parts by writing 


1(8)=1q(9) + 1,(9) 


where J,(@) is the scattering function due to diffraction and J,(0) is that due to 
reflection and refraction. The former predominates at small elongations and 
describes the F-corona, and the latter at large elongations and describes the 
zodiacal light. For the first, Fraunhofer diffraction theory leads to the result 
J ,?(x sin 6) 

sin? 6 (3) 
where x=ka=2za/A. On the other hand, the variation of the brightness of the 
zodiacal light with elongation does not depend strongly on the form of the function 
1,(@), and for convenience we write 


71,(8)= ty (4) 


which is the form for an isotropic scatterer with albedo y. That there is some 
physical justification for this simplification is shown by Richter (5) who examined 
the scattering by particles of various sizes in the laboratory and found that for 
dielectric particles of diameter less than 10~* cm the scattered intensity is approxi- 
mately independent of @ for values of 6 greater than about 40°. 

We have considered two.different functions n(a). The first, which has been 
used by van de Hulst (3) and Beard (6) is 


n(a)da=Ca~? da (5) 


nl (6) = 


and the second has the form 
n(a)da= Dae~'*" da (6) 

where C, D, p and aare constants. In the first distribution the smallest particles 
are the commonest, whereas in the second the most numerous particles are 
distributed over a range of sizes rather than being concentrated at the small end 
of the distribution. The first function has the merit of simplicity and appears to 
describe the distribution of micrometeorites (see Section 7(ii)), whereas the second 
might describe more accurately a mixture of distributions, e.g. of particles of 
different materials, each with a different value of a. 

Considering first the diffracted part of the scattered light and the first of the 
above functions, we have, from equations (2), (3) and (5), and since diffraction 
occurs in the forward direction only 


J ,?(xsin @) 
(x sin 6)?-* 


a, fn/* 
F (cs k)=J,dAR cosec!**e | | Cr? dé da. 


J ,*(kasin @) 
(ka sin 6)?-* 


n/* (as 
=J,d\R cosec!*e | | j2-aqi-a—p dado. 











160 M. F. Ingham 








Putting ka sin 0 =2, da= id , we have, after some reduction, 
ksin@ 
ar |? 
Fa (e ;k)=J ,dXARCk?— cosec!** € | sin*t?—>6.G,,(0; k) dé (7) 
where 
: ka, sin Lm 2(z) ‘ 
or 
Fy, (€;k)=J dARCR” * cosec!**. H, » (€;R) (8) 
where 
/* 
H, ,(€;k)= | sin*+?—59 , Gp(0; k) dé. (8a) 


For the light scattered by reflection and refraction we have, using equation 
(4), and because J,(@) is independent of a, 


F, (¢)=J, dAR cosec!+* | "by sin® 0 d0 | “Cat? da (9) 
or, for p¥ 3, | 
F, (c) =u ,dARCy cosec!+*e . H,,(e) (se) (10) 
where 
H,,(e)= { "} sin 640. (104) 


Now a, >a, and (Section 7(ii)) on the basis of observations of small meteors, p 
will be given the values 4 and 5. In this case we may ignore a,°~” in comparison 
with a,°~” and in place of equation (10) write 





3—p 

F, (€) =J ,daARCy cosee!%e. H,, (6) (A -). (11) 
The mathematical details of the steps resulting from the use of equation (6) 
are similar to those just given and since what follows applies equally to both forms 

of n(a) we shall consider only the first, giving the final results for both at the end. 
In order to discuss the polarization of the light scattered by small particles it is 
necessary to know, or assume, something about their nature. The only certain 
samples of interplanetary material to hand are meteorites, and the work of 
Kloverstrom and Rense (7) on the polarization of light reflected from the plane 
ground surfaces of meteoritic samples, and also the work of Richter already referred 
to (5), indicates that the degree of polarization, p(@), varies with the angle of 

scattering, 0, as a function of the form 
sin? 6 

p(@)= q—sin®0 (12) 
where q is a constant depending on the nature of the sample, for example whether 
it is stony or metallic. 
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If we adopt for p(@) a function of the form given by equation (12) then we 
may use equation (g) to write the expression for the polarized component of the 
dust scattered light as 


f,,(€) =J,aAR cosec!**e | “byp(6) sin=6 d0 | ” Ca?-eda 
€ a 


sin2+# @ 


=J,dn peg gy SE 
,4AR cosec I ere 


as 
d8 | Ca??? da. 
ai 
Thus the observed polarization of this component is given by 
* sin***@ 
) me [>a 
te 


f sin* 6 d@ 


€ 


Pr)= 2 





(13) 


4. Criteria for the choice of the electron density distribution.—Once a model of the 
electron density has been chosen, K and p, can be computed at all elongations from 
the equations given by van de Hulst (2) and, from the observed values of py and 
K+F, values of F and p, are derived. Lastly, F may be further separated into 
a component due to diffraction, F,, and one due to reflection and refraction, F.. 
By analogy with the total corona we have, in an obvious notation, 


pp F=p,(FatF,) =Pr Fa + py,F,. 
Here, though, py, 18 zero since we have used Fraunhofer diffraction theory in 


deriving an expression for Fy, and so 


p= Be 
Pr, 

A theoretical expression for p,, was obtained in Section 3 (equation (13)) and 
was seen to depend both on the value of « and on the nature of the scattering 
particles. However, for the range of values of « and g with which we shall be 
concerned, p, is constant for « < 10° and so an error in assigning a value to p,,. 
does not affect the manner in which F, varies with elongation. Finally by sub- 
traction of F’, from F we obtain Fy. 

Two criteria for the choice of n,(r) are available. The first concerns the 
variation of F, and F, withe. From equations (8) and (10) we see that 

Fy(e)xcosec'**e. H, ,(€) (14) 
and 

F,(€) ccosec!*e. H, (€). (15) 
Now H,,,(e) and H, (€) vary only little with « as compared with cosec!+*e, The 
slopes of the graphs of log F,(¢) and log F,(e) as functions of log « are therefore a 
measure of the value of « in the regions from which come the major contributions 
to the functions H, ,, (€;k) and H, (€), which, it will be remembered, are integrals 
evaluated along the line of sight. 

Here it is necessary to anticipate the contents of Sections 6 and 7 and to say that 
two values of « will be considered, «= 1 and «= 3/2 and two values of p, p= 4 and 
p=5. Consider for example the case « = 3/2, p= 5 and n(a) given by equation (5). 
For a wavelength of A= 6283A (k=27/A= 10°) the function sin*+”—> 6. Gp(0; k), 











162 M. F. Ingham Vol. 122 


whose integral between the limits « and 7/2 yields the function H, ,(€;), has 
such a form that for «= 1°, say, half the integral comes from the region defined by 
1° <6< 14° and half from the region 14°<@<go°. Furthermore, 10 per cent of 
the integral is contributed by the region 1° << 4° and 10 per cent by the region 
30°<@<go°. Inan acceptable sense, therefore, 80 per cent of the integral comes 
from the region defined by 4°<@<30°. Slightly different limits hold for other 
values of e«, and Table I shows the distances, in solar radii, which these limits are 
from the Sun. The function }sin*0, whose integral from « to z yields the func- 
tion H, (€) is symmetrical about 6 = go° and 80 per cent of the integral comes from 
the region bounded by @= 40° and 6= 140°, whose limits are situated at distances 
from the Sun which are also given in Table I. 


TABLE I 


Elongation, scattering angle and distance from Sun. The 
values of 0 define the range from which comes 80 per cent 
of the functions H3,2, s(€; k), k= 10° (a), and H,,, (€) (b). 





€ 6° rin Ro 
(a) I 4 54 
3° TS 
3 5 86 
30 15 
3 5 128 
30 22 
4 6 143 
30 30 
5 7 153 
30 37 
(d) I 40, 140 6 
2 40, 140 12 
3 40, 140 18 
4 40, 140 23 
5 40, 140 29 


Thus the form of these functions is such that for any given elongation <5° 
the major part of the component F, arises from scattering by particles which lie 
nearer to the Sun than those which give rise to the component Fy. If therefore 
we find, for a given n,(r), that log F’, varies linearly with log « as far as «= 5°, then 
we can deduce that/out to about 0-5 A.U. « has the constant value given by the 
gradient of the line, and hence that F,, which arises closer to the Sun, must vary 
with ¢ in the same way as does Fy. If it does not, i.e. if the graph of log F, against 
log « either is not a straight line or has not the same slope as that of log F,, then 
the model chosen for ,(r) is not satisfactory. 

The second criterion concerns the function p,. At elongations greater than 
40° F, is small in comparison with F,andp,~py,,. With decreasing « the rapidly 
increasing effect of the unpolarized component F, results in a sharp decrease of p,» 
relative to p,. At very small elongations (« <2°), where the functions H,, ,(¢) 
and H, («) are approximately constant as compared to the behaviour of cosec'**e, 
equations (14) and (15) show that the ratio of F, to F, is very nearly independent 
ofe. Itis therefore to be expected that, with decreasing ¢, p » will decrease rapidly 
at first and finally approach a constant value close tothe Sun. At any rate, there 
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is no reason to suppose that p, will increase with decreasing ¢ at any point, and 
hence a model for n,(r) which indicates that it does cannot be regarded as accept- 
able. Ohman (8), at the eclipse of 1945 July 9, measured the polarization of the 
F-corona and concluded that it was very small in comparison with that of the 
K-corona. In fact he found a small negative value at an elongation of o-19R_, 
but admitted that the material on which his analysis was based was ‘‘ very meagre ’’. 
The strong positive polarization of the zodiacal light makes it improbable that the 
F-corona should be negatively polarized. 

In this argument we have ignored the fact that there must be a zone free of dust 
surrounding the Sun within which all particles are vaporized. Such a zone would 
produce a more or less sharp cut-off of F,, while affecting F, to a much lesser 
extent, and this in turn would cause a cut-off in p,. Over (g) considers that a 


probable value for the radius of this zone is 4R, and its presence in no way 
affects the argument that p, must decrease with decreasing «. 

5. Construction of an observational model of the corona and the zodiacal light.— 
In order to apply this method of separation we require a set of measurements of 
K+ and py both for the corona and the zodiacal light, and here we are confronted 
by a difficulty for no one observer has made observations of the surface brightness 
and polarization of the corona from its innermost to its outermost parts. 
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Fic. 2.—The “ observational’’ model of the corona and zodiacal light-surface brightness. 








164 M. F. Ingham Vol. 122 


Van de Hulst (2) has, on the basis of earlier observations and on the assumption 
that p,.=0, constructed a model corona, both for sunspot maximum and for 
sunspot minimum, as far out as x=6R_. In the outer corona there are the 
observations of Blackwell (10) who measured K + F in the range 6 <x < 55R_ and 
Po in the range 6<x<20R,. In region of the zodiacal light the most recent 
observations are those obtained at Chacaltaya which extend over the range of 
elongations 20° <«<70°. 

The difficulty alluded to above arises from the fact that the van de Hulst 
models do not join on smoothly to Blackwell’s observations, though these can 
reasonably be extrapolated to link up with the zodiacal light measurements (see 
Fig. 7 of Paper I in this series). In constructing an observational model of the 
corona and the zodiacal light we have therefore reconciled the inner two sets of 
observations by using the values of the brightness in the van de Hulst maximum 
model out to an elongation of 2K, and then interpolating so as to join on smoothly 
to Blackwell’s observations. It is true that Blackwell’s observations were made 
at sunspot minimum (1954), but the van de Hulst minimum model is even less in 
agreement with them than is his maximum model. We have altered the outer 
values in the van de Hulst model rather than the inner values in Blackwell’s 
range of observation because we consider the latter to be the more reliable since 
Blackwell observed from an aircraft flying at 30000 feet and had vastly better 
sky conditions than are obtainable at a ground station. 

Fig. 2 shows the resulting ‘‘ observational ’’ model of the surface brightness, 
K+F, of the corona and the zodiacal light at a wavelength of 6200A. The 
continuous curve represents the model. The triangles are van de Hulst’s 
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Fic. 3.—The ‘‘ observational’’ model of the corona and zodiacal light-polarization. 
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ion maximum model, the open circles are Blackwell’s coronal observations and the 
for filled circles are an average of those for the zodiacal light. 
the Fig. 3 gives the values of the polarization, pp, in the observational model. 
and From x= 1R_ to x=2R,,, the range in which we are prepared to allow the assump- 
ent tion p,=0, these are deduced from the equation 
> of 
Po= - K 
ulst A+ 
can and the values for the quantities given by van de Hulst. These points are marked 
(see by the triangles. The open circles are again Blackwell’s coronal observations 
the and the filled circles those for the zodiacal light. Interpolation bridges the gaps. 
s of 6. Separation of the model.—The choice of a function n,(r) which should 
um satisfy both the conditions laid down in Section 3 proved to be a delicate matter 
thly since the variation of both F, and p, with « depended strongly on the model 
ade chosen. The form of ,(r) selected was 
ss in Seine 
ater n,(r)= ge (r>6R_) 
ell’s ’ 
ince and Table II gives the values of n,, Fg, F, and p, for this model. The value of « 
‘tter indicated by the gradients of log F, and log F, is unity, or, if the values at smaller 
ess, Tane II 
The 
Ist’s Values of n,(x), F 4, F,, and pp for the selected model of n,(r) 
x Fg F, 
i RO n,(x) 10-!4BO 10-!4BO Dr 
6 31300 736 42°3 0°016 
8 13200 362 22°0 0°017 
10 6750 211 13°5 0-018 
12 3910 144 9°83 0°019 
14 2460 106 7°22 0-019 
16 1650 80°1 5°64 o'o1g 
18 1160 63°6 4°71 0°020 
20 843 52°38 4°09 0°021 
— elongations are considered to be the more reliable, « = 3/2 is admissible. Linking 
this function smoothly to the values of the van de Hulst maximum model in the 
range 1<r<4R_ completes the model of the electron density. Separation of 
K+ F into K and F can now be effected at all elongations and values of p, and p, 
derived. Columns 3 to 8 of Table III show the results of the first separation. 
At elongations greater than 55° H,, ,(€;&) (k= 10°) vanishes for «= 1 or « = 3/2 
and p=4 or p=5. We then have the situation in which K is negligible in com- 
parison with F and F, is negligible compared to F,. This means that for e< 55°, 
Po=Py=Py, Thus observations of the zodiacal light at large elongations are 
— observations of F,, and these measurements (Fig. 2) lie nearly on a straight line of 
sine: gradient —2-4. Inthe range 20° <«< 70° the graph of log {cosec®e . H,,,(€)} asa 
function of log « approximates to a straight line of gradient — 2:4 and so we can 
deduce that, for r> 0-5 A.U., «=3/2, at least to 1 A.U. 
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TABLE 


Separation of model 





Elongation 





x/RO € 


K+F 
10-!4BO 


K 
10-!®?BO 


F 
10-!2BO 


PK 





Fare) 16’ 


I'l 
I'2 
1°4 
1°6 
1°8 
20 


5°54 x 10° 
1°74 X 108 
7°36 x 10° 
2°20 x 10° 
9:06 x 104 
4°47 x 108 
2°52 Xx 104 
8-61 x 10° 
5°43 x 10° 
2°40 X 10% 
1°35 x 10° 
851 
407 
234 
159 
115 
87°1 
69°2 
57°5 
37°2 
25°1 
17°8 
12°6 
9°33 
6-92 
5°25 


2°35 

1°38 

0884 
0*609 
0°442 
0°332 
0°258 
6-204 
0°166 
0°137 
O°1I5 


5°34 x 10° 

1°62 x 10° 
6°62 x 10° 

1°84 x 10° 
6°97 x 104 
3°09 X 104 
1°52 x 104 
3°31 x 108 
1°55 X 10° 

392 
150 

72°5 
22°9 

9°38 

4°53 

2°44 

1°43 

0°89 

0°59 

0°24 

O'12 
0°062 
0°037 
0°023 
O'OI5 

oo) fe) 


0°003 
0-oo1 
000! 
0°000 


1°99 x 10° 
1°16 x 10° 
7°40 X 104 
3°56 x 104 
2°09 X 104 
1°38 x 104 
1°00 X 104 
5°30 x 10° 
3°88 x 10 
2°01 x 108 
1°20 X 10 
778 
384 
225 
154 
113 
85°7 
68-3 
569 
37°0 
25°0 
17°7 
12°6 
9°31 
6-90 
5°24 


2°35 

1°38 

0884 
0‘609 
0°442 
0°332 
0258 
0°204 
0166 
0°137 
O'1I5 
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We can now compute p, from equation (13), choosing q so that py =pPp 


Table IV shows the results of this calculation for «=1, g=3°2 and 
a=3/2, g=3°3. For «>20° we require the values in the last column while for 
x <20R,, we have seen that there is a case for either value of «. F, may now be 
obtained by dividing p,F by py, Fy by subtraction of F, from F and the separation 


is complete. The functions p(0) which result from these choices of g are tabulated 


at «e=60°. 
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corona and zodiacal light 





Pr, F,, 10-"*BO Fa, 10-*BO 





a=1 a= 3/2 a=I a= 3/2 cm~* 





1°99 x 10° 1°99 x 10° 4°03 x 108 
1°16 x 10° 1°16 x 10° 1°60 x 108 
7°40 X 104 7°40 X 10 7°08 X 107 
3°56 x 104 3°56 x 104 2'27 X10" 
2°09 x 104 2°09 X 104 I‘OI X 107 
1°38 x 10¢ 1°38 x 104 5°19 X 10° 
1°00 x 104 1‘00 x 104 2°81 x 10° 
5°26 x 10° 5°27 x 10° 7°59 x 108 
3°81 x 10° 3°81 x 10° 3°98 x 10° 
1°93 X 10° 1°93 X 10° 1°17 x 10° 
1°14 X 10° 1°14.X 10° 5°50 x 10 
734 735 3°13 x 108 
361 361 1°32 x 104 
210 211 6°75 x 10° 
143 144 3°91 x 10° 
105 106 2°46 x 10° 
79°9 80°1 1°65 x 10 
63°4 63°6 1°16 x 10% 
52°6 52°8 843 
34°1 34°2 433 
22°8 22°8 250 
15°9 15°9 
II‘I II‘! 
7°92 7°97 
5°68 5°72 
4°15 4°19 


1°56 

0680 
0°247 
0°073 
0'025 
0006 
0002 
0'00I 
0:000 





























in Table V. Columns 11 to 14 of Table III exhibit the results of the second 
separation. 

Fig. 4 shows po, Px, Py and py, the last for «=3/2. The relations of the 
quantities may be understood as follows. At large elongations K+ F= F =F, and 
bo=Pyp=Py,- Asweapproach the Sun weat first still have K+ F= Fandp,)=p,, but 
py diminishes in relation to pg, because of increasing dilution by the unpolarized F4. 
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TaBLe IV 
Values of p F, 


¢ Pr, 
a=1 a= 

° 0'279 0°290 
10 281 ‘290 
20 287 "295 
30 296 *300 
40 “306 ‘309 
50° “315 *317 
60 *320 "324 
7° ‘317 324 
80 "303 312 
go 279 ‘290 
100 246 °255 
110 ‘203 213 
120 "159 ‘169 
130 115 123 
140 076 “082 
150 "045 O51 
160 018 "029 
170 "007 003 
180 0000 0°000 


Closer to the Sun at about 60R,, we begin to have py>py, as the strongly 
polarized K-component becomes effective, and at about x=25R_, py begins to 
increase again and approaches p,,, though we never have py =p, because of dilution 
by F;. Meanwhile p, at first approaches a constant value, in sympathy with the 
behaviour of p,,, and then decreases sharply when the line of sight begins to pass 
through the dust-free zone. 


TABLE V 
Velue 6 te in eh 0 ae 
alues un on aetemganen i, 
ati q—sin? 0 9 
e° p(9) 6° 
a=1, q=3'2 a= 3/2, 9=3°3 
° 0*000 0°000 180 
10 ‘ 009 “009 170 
20 038 "037 160 
30 “085 082 150 
40 148 "143 140 
50 "225 216 130 
60 -306 "294 120 
70 381 +365 110 
80 "435 “416 100 
go 0°455 0°435 go 


7. Albedo, mass density and number density of the solid particles.—Betfore 
we can evaluate the expressions in equations (8) and (11), developed in Section 3, 
we must consider what values to assign to the constants a@,, a2, p and a); values of 
a have already been decided upon. We can then use the ‘‘ observed ”’ values of 
F, and F, to deduce C and y. 
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(i) The value of a,. The lower limit to the size of particle present in inter- 
planetary space is set by the point at which the repulsive force of radiation pressure 
balances the gravitational attraction of the Sun. This limit depends on the 
density of the particles and if it be allowed that the material composing the 
zodiacal cloud is similar to that of which meteorites are made we may estimate this 
quantity. ‘The majority (> go per cent) of meteorites seen falling are stones (11) 
with a density of about 3gmcm~*. On the other hand, Whipple (12) suggests 
p=2 for an aggregate of stony-iron particles. Moreover, the finest dust may 
possibly originate from porous friable material of low density. We shall therefore 
consider the two cases p = 2 and p = 3 as being representative figures for the density 
of the smallest particles composing the zodiacal cloud. 

For black spheres the lower limits are a, = 0-3 and a, = 0-2, in the two cases. 
The problem of the pressure of radiation on small, totally reflecting spheres has 
been discussed by Proudman (13) and use of his results leads to the values a, = 0-4 
if p= 2 and a, =0°3p if p= 3 at a wavelength given by k = 27/A= 10°, i.e. A= 6283 A. 
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Fic. 4.—Polarization of the components of the “‘ observational’’ model. 


(ii) The value of p.—The observations of both Watson (14) and of Lovell 
(15) on micrometeorites yield a value for p of about 5, though these results refer 
to particles whose radii lie in the range 100% <a< 1000u, whereas the present 
investigation considers particles whose radii are 11 and less. Since other authors 
(3, 6) have found smaller values for p, we have considered the two cases p= 5 and 
p=4. 

(iii) The value of a,. ‘This is very uncertain but must be several orders of 
magnitude larger than a, if the distribution defined by (5) or (6) holds for values 
of aup to1mm (1000). In Section 3 we saw that we may treat the upper limits 
involving a, as infinite for the values of p selected above. 
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The integrals in equations (7a) and (8a) were evaluated graphically at a wave- 
length of 6283A so that the expressions for Fz and F, could be compared directly 
with the observations of the outer corona and the zodiacal light made at an effective 
wavelength of 6200A. The functions 


log {cosec'**e. H, ,(€; k)} and log {cosec!**e. A, (€)} 


tonit2 +4 Ro 


log cosec’"€.H_ (ek) 











T T T T ee T a 
4 30 40 506070 
€ in degrees 


Fic. 5.—The function cosec 1+%e, Hy,» (ek) for p=5, a,=0'3n, R=10° and a=1 
P and a=3/2. 


were then plotted against log « on the same diagrams respectively as log F', and 
log F,, (Figs. 5 and 6) and the values of 


3-—p 
(J,dARCk?-*) and (J. darcy a4 -) 


deduced by securing the best agreement between the computed functions and the 
observations. Values of C and y were then obtained. 

The mass density and number density of particles in interplanetary space neat 
the plane of the ecliptic may now be derived. Using equation (5), the total mass 
in grams, M, of the dust particles in 1 cm* at 1 A.U. from the Sun is given by 


M= tmp | a°n(a)da 


a, 
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and the total number N, of particles of all sizes in 1 cm* at 1 A.U. is given by 


N= a n(a) da. 


Since C is known the evaluation of these expressions is elementary. 








log,.cosec’e.H, €) 
log,,cosec"e.H,(€) 
"a 


109, 








UJ ‘a 


T —~ «tf T ae 
4 6 20 30 40 506070 
€ in degrees 


. 6.—The function cosec 1+%e. Hy,(¢) for a=1 and a=3/2. 


TaBLe VI 


Albedo, mass density and number density of interplanetary dust particles. 
M and N refer to the plane of the ecliptic at 1 A.U. from the Sun 


First distribution: n(a) da=Ca~? da 
ay Y c M 

10-4 gm cm~* 

a=1 a=3/2 a=1 a= 3/2 @2=1 a= 3/2 

0719 048 3°9X107%* 1:6x107%* = = 5-1 2'0 

0°23 0°56)=—_ 20 X 10-8 7-8x10731 0-82 0°33 

0°16 0°99 «66 -a x 10-**' agxzo—™ §¢:3 a1 

o18 os1 4°3X1078° 1°6x107%9 ogo 0°33 


Second distribution: n(a) da=Dae~®/%* da 
A Qa, 9 4 D M 
10cm 10-*cm 10~*4 gm cm~% 
a=1 a=3/2 a=1 a=3/2 a=1 a=3/2 
0°3 0716 ov51 grx ro. gx t0"* a3 (0-71 
03 0°16 = 0°35 179x107 87x10" 48 2°2 
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The results of these calculations in the various cases considered are shown in 
Table VI. 

No account has so far been taken of the effects of the finite size of the Sun or of 
the dust-free zone surrounding it. The former will cause an increase in the 
values of F’, over those calculated assuming a point source. ‘The latter will cause 
F.. to decrease sharply as soon as the line of sight begins to pass through the zone; 
F, will be less affected since it arises from particles more distant from the Sun 
than those responsible for F,. 

The first effect was treated by the method described by van de Hulst (3) and 
the second column of Tables VII and VIII shows the amount by which the value 
of log,» /4(x) calculated from equation (8) must be increased in order to take 
account of it. 


TasB_Le VII 


Change in log,) F, due to the finite size of the Sun and the presence of a dust-free zone 
of radius 4Ro—first distribution 


Increase 
due to 

finite Decrease due to dust-free zone Total change 
size of 

Sun 


** Observed” 


increase 





allp,a, 274 P=4 P=5 P=5 P=4 P=4 P=5_ p=S 
*) Gy =0'3p A, =0'4p A, =0'3 A, =O0'4p A, =0'3p A, =O0'4p a, =0'3m a,=O'4p 





0'287 0184 0460 # 0°286 
0°326 ‘246 “158 "404 ‘243 ? +0168 —0'078 -+0°083 
‘232 "219 136 *359 ‘208 + °* + 096 — ‘127 + ‘024 
-146 "174 ‘097 "284 "144 ‘ + 049 — ‘138 + ‘002 
"103 139 ‘070 ‘227 ‘101 ‘ 033 — ‘124 + ‘002 
078 “E33 "O50 178 068 -—°: t °o28 — ‘100 + ‘o10 
‘061 ‘087 "035 +138 "O45 : 026 — ‘077 + ‘016 
034 ‘037 ‘009 064 ‘009 ‘ "025 "030 -++ ‘025 
025 018 ‘000 035 70000 ++ * 025 ‘O10 + ‘025 
‘O15 ‘000 ‘000 - 4 r O85 + OLS + “Ors 
‘009 ee + ‘009g "009 + ‘009 
‘006 - 9 + ‘006 + ‘006 + ‘006 
‘003 a, r @3 + OCs + O03 
‘000 + ° + ‘000 + ‘000 + ‘000 


I’ 
I 

I 

I° 
1° 
1° 
2° 
2: 


> wW 
° 


*165 
‘070 
048 
"036 
032 
026 
“000 


Recently Over (g) has considered theoretically the vaporization of solid 
particles of SiO, near the Sun and concludes that half the particles would have 
evaporated at 4 solar radii from the Sun’s centre. Here we shall consider a crude 
model in which there is a sharp cut-off in the density of solid particles at 4R>. 





served” 
crease 
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There will thus be a cut off in the integral which yields H, ,(€; &) and instead 
of equation (8a) we must, in the range 1<x<4R_ write 


% 
H,, ,'(;k)= | sin*+?-5 0. Gp(0; «)d0 


where 0, is the scattering angle at the point where the line of sight enters the dust- 
free zone. ‘lables VII and VIII show the amounts by which log, /4(x) must be 


TasLe VIII 


Change in logy F 4 due to the finite size of the Sun and the presence of a dust-free zone 
of radius 4R~>—second distribution 


Increase due 
x to finite size Decrease due Total change ‘‘ Observed” 
Ro of Sun to dust-free zone increase 


all a A=1 ph Ag=2h aAg=1p ag=2p 


0082 = 0022 
0°326 ‘O70 ‘017. +0°256 +0°309 
°232 "064 ‘O12 -168 + °220 
*146 048 ‘006 098 "140 
"103 037 ‘003 066 *100 
‘078 "029 ‘OO! "049 ‘O77 
‘061 "022 ‘000 "039 ‘061 
034 ‘O10 "024 034 
025 "004 ‘O21 025 
‘O15 “000 ‘O15 ‘O15 
‘009 ‘009 ‘009 
‘006 ‘006 ‘006 
003 003 003 
“000 *000 *000 


t+ettetete+ett¢et 
t++tetet+++++ 
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decreased in consequence in the various cases and also gives the total change in 
log,» /,,(x), while the last columns show the increase in the ‘‘ observed ’’ values of 
log F', over a straight line of slope —(1+«) passing through the point which 
represents log F (1°) (Fig. 5). 

8. Discussion of results.—Although the results of the preceding section are 
based on a smoothed mean of a large number of different observations it is worth 
while discussing the results obtained from the model. ‘There we will consider 
under five headings. 
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(i) The values of y, M and N.—The values of the albedo y in ‘Table VI are 
consistently higher for «=3/2 and for «=1. Kloverstrom and Rense (7) 
obtained values of 0:27, 0°62 and 0-38 for three stony meteorites and 0-58 for a 
metallic sample. For either value of « it is clearly not possible to select either 
the distribution of particle sizes, or a value for p, a, or ay on the basis of the value 
of y. 

The mass densities are of the same order of magnitude as those deduced by 
Behr and Siedentopf (16). Recently attempts have been made to measure the 
density of solid particles near the Earth with satellites and space probes, but 
although the results appear to agree with ours to within one or two orders of 
magnitude, they differ widely among themselves. Since no other independent 
measurements of this quantity exist the differences between the values of M in 
Table VI afford no criterion for selecting values of p, a, or dp. 

The values of N are from one to three orders of magnitude greater than 
previous estimates, the reason being that we have accounted for the zodiacal light 
and the F-corona in terms of many small particles rather than fewer larger ones. 
Behr and Siedentopf (16), Blackwell (17) and Allen (18) considered models 
involving particles of one size only and those sizes were one to three orders 
of magnitude greater than the values of a considered here. 

The values of both M and N may be compared with the figures for the inter- 
stellar grains near the galactic plane (19) which are M=1-4 x 10-6 gmcm~ and 
N=200km-*, These values are, however, derived on the assumption that 
p=1'tgmcm~* and a=3 x 10-°cm. 

(ii) Surface brightness of F-corona close to the Sun.—Tables VII and VIII 
show that only the first distribution with «=3/2, p=4 and a=o-4yp and the 
second distribution with « = 3/2 and both values of a) account at all satisfactorily 
for the rise in log F', close to the Sun; the results for «=1 are in all cases bad. 
The figures for the ‘‘ observed ’’ increase given in the last columns are, however 
somewhat arbitrary. 

(iii) Number of particles in the range of sizes 10-* <a < 10-1! cm.—Van de Hulst 
(3), who considered the distribution of particle sizes given by equation (5), arrived 
at the values C= 10~*°, p= 2-6 and comparison with Watson’s analysis of meteor 
frequencies (14) led to the conclusion that ‘‘the number of particles with sizes 
between 1mm and o-1 mm in interplanetary space is about 10000 times larger than 
can be inferred from the number of telescopic meteors’’. Using the values C= 10, 
p =2°6, a,=10~*cm and a,= 10~' cm the expression 


("le 
| Ca” da 


Ay 


has the value 10~!’, whereas our values of C~ 10-*, p= 5 and C~ 10~*6, p = 4 lead 
to values for the integral of the order 10~** and 10~*° respectively. The discre- 
pancy between the predicted number of telescopic meteors and those actually 
observed is thus to a large extent removed. 

(iv) The polarization function p(@).—'Table V shows that in order to account 
for the observed polarization of the zodiacal light in terms of scattering by dust 
particles alone we have had to assume a function for the degree of polarization 
due to a single particle which has a maximum value of about 45 per cent when the 
angle of scattering isgo”. Sucha large value has never been obtained in laboratory 
measurements. Kloverstrom and Rense (7) obtained maxima of about 30 per cent 
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for stony and 20 per cent for iron meteorites, though this work dealt with large 
plane surfaces and not with a cloud of tiny particles. Richter (5) worked with an 
aerosol of iron particles prepared from iron carbonyl and also with a suspension of 
small particles of quartz in a sealed room and recorded maximum polarizations of 
15 per cent and 30 per cent for dielectric and metallic particles respectively. It 
is, however, doubtful just how far the materials and physical conditions in these 
experiments are representative of interplanetary space. 

It is possible that independent evidence may be obtained from the study of 
noctilucent clouds. ‘These are observed mainly in high northern latitudes and 
are formed at an altitude of between 80 and go km, but the nature and origin of the 
particles of which they are composed is uncertain. Witt (20) has measured the 
polarization of the sunlight scattered by such clouds. ‘The measurements were 
made in red and blue light and at scattering angles of up to 60°. The degree of 
polarization exceeds 40 per cent at 6 = 60° in both colours and appears to reach even 
higher values at larger angles of scattering. He interprets these measurements 
according to the Mie theory in terms of scattering by dielectric spheres, and shows 
that the curves indicate a particle radius of 0-1—o-2 4, depending on the refractive 
index. Without asserting any connection between these particles and the inter- 
planetary dust we nevertheless have here a mechanism capable of producing an 
even greater degree of polarization than appears to be required to account for the 
observed polarization of the zodiacal light. 

(v) The reddening of the corona.—We consider finally the variation of Fy 
with k(=27/A). It was predicted by van de Hulst (3), observed by Allen (18) 
and confirmed by Blackwell (17) that the F-corona is redder than the Sun. 
Van de Hulst predicted that F, should vary with k as k-** and Blackwell showed 
that the reddening of the total corona increases with elongation as the reddened 
F-corona increases relative to the almost colourless K-corona. 


TaBLe IX 
The reddening of the corona 
1st distribution 2nd distribution 
p=4 p=5 4=104%cm ad=2xX107*cm 


a=I 1°02 1‘00 0'99 1'02 
a= 3/2 1°13 1°13 1°09 1°12 


This reddening may be described by the ratio 


hui I (Ag, %2) | L(Ag *1) 
I(Ay ¥2) | L(Ay #1) 





where /(A,x) is the surface brightness of the total corona at wavelength A and 
elongation x. Blackwell, using the values A, = 4300 A, Ag=19 000A; x,=1°5R_, 
%, = 2°5R. obtained R= 2:17. It would not be correct, though to use the analysis 
developed here to calculate a theoretical value for this ratio since A= 2-9 x 10-4 cm 
and a,=03 x 10-*cm leads to x,=27a,/A~1, for which Fraunhofer diffraction 
theory is certainly inapplicable. Allen, however, measured this ratio at wave- 
lengths 4000 A and 6500A and at elongations 12K, and 2°4R,, and obtained 
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the value 1°48, though Blackwell considers this to be too high. On the other 
hand Ney, in a private communication, states that at the eclipse of 1959 October 2 T 
he obtained the value R = 1-08 in the equatorial region for A, = 4750 A, A, = 8300 A; 
%=12R_, x2=2:1R... Table IX shows the values for a,=0-4x 10cm; 
A= 4189A (ky=1'5x 10), Ay=6283A (kg=10°); xy=1'2RO, %2=2-4R_. 
There is a clear distinction between the values for «= 1 and those for «= 3/2 and 
the latter are in reasonable agreement with Ney’s result, though considerably 
lower than Allen’s. Further observations of this quantity are required. 

Our treatment thus appears able to account for the observed trends in the 
F-corona and the zodiacal light and the figures given in Table VI are plausible. 

g. Acknowledgments.—My thanks are due to Professor D. E. Blackwell and 
Professor R. O. Redman for much advice and many fruitful discussions of this 
problem. I am particularly grateful to Professor E. P. Ney for permission to 
refer to a result of his work prior to its publication in the Astrophysical Journal. 
I am indebted to the Electors to the Isaac Newton Studentship for support 
during the period in which this work was done. 

The Observatories, 

Madingley Road, 
Cambridge: 
1960 October to. 


( 

References L 
(1) D. E. Blackwell and M. F. Ingham, M.N., 122, 129, 1961. 
(2) H. C. van de Hulst, B.A.N., 11, 135, 1950. v 
(3) H. C. van de Hulst, Ap. F., 105, 471, 1947. 
(4) H. C. van de Hulst, Light Scattering by Small Particles, Chapter 8, Wiley, 1957. (< 
(5) N. Richter, Mon. Deut. Akad. Wiss. Berlin, 1, 727, 1959. 
(6) D. B. Beard, Ap. 7., 129, 496, 1959. ( 
(7) F. A. Kloverstrom and W. A. Rense, Ap. 7., 115, 495, 1952. ( 


(8) Y. Ohman, Stockholm Obs. Ann., 15, No. 2, 1947. 

(9) J. Over, Proc. Kon. Ned. Akad. van Wetenschappen, 51B, 74, 1958 ( 
(10) D. E. Blackwell, W.N., 115, 629, 1955. ti 
(11) F. Watson, Between the Planets, Harvard Books on Astronomy, 1941. 

(12) F. L. Whipple, Ap. 7., 111, 375, 1950. ( 
(13) J. Proudman, M..N., 73, 535, 1913. e1 
(14) F. Watson, Harvard Annals, 105, 623, 1937. ( 
(15) A. C. B. Lovell, Meteor Astronomy, p. 139, Oxford, 1954. 
(16) A. Behr and H. Siedentopf, Zeits. f. Astrophys., 32, 19, 1953. 
(17) D. E. Blackwell, M.N., 112, 652, 1952. 

(18) C. W. Allen, M.N., 106, 137, 1946. (¢ 
(19) C. W. Allen, p. 226, Astrophysical Quantities, London, Athlone Press, 1955. 
(20) G. Witt, 7. Geophys. Res., 65, 925, 1960. 








oA; 
cm; 


and 
‘ably 


. the 
ible. 

and 

this 
n to 
rnal, 
port 


57: 











THE K-CORRECTION AND “‘STEBBINS-WHITFORD”’ CORRECTION 
IN RELATIVISTIC COSMOLOGY 


R. Van der Borght 
(Communicated by G. C. McVittie) 


(Received 1960 November 2)* 


Summary 


In the present paper we derive an expression for the K-correction and 
the ‘‘ Stebbins-Whitford ’’ correction in relativistic cosmology, correct to 
the second order in the redshift 5, and with coefficients calculated for the 
present epoch. The results are then compared with those derived by 
McVittie and Davidson and it is shown that, if one works to this order of 
approximation, the two corrections cannot be completely separated and 
that a “‘ mixed ”’ correction has to be introduced. 





In the study of a uniform model universe, one is led to the following formula 
(McVittie 1956, p. 156, Davidson 1959, p. 60) connecting the luminosity distance 
Dand the apparent magnitude m of a given galaxy 

log) D=0'2 (m— M,)+1+F (1) 


where 


(a) F=} log, scar | : o(A)B (i 4 - 5) a (2) 


(b) M, is the absolute magnitude of a typical galaxy 

(c) 4, is the ‘‘ Hubble constant ”’ 

(d) +=t)—t is the time of travel of light emitted at the time ¢ and observed at the 
time fy. 

(e) B(h,7,A) dA is the energy emitted by a galaxy at the time ¢ in the range dA in 
ergs/sec cm”, 

(f) (A) is an empirical function depending on the atmospheric extinction and 
the response of the apparatus. 


(e) i= | “ o(A) B(o,A) da. (3) 


Writing 





A 
—_} = —A(6-& 
B(x : 5) Bix, A—A(8—82)} 
where 
x=h,t 
we have, using a Taylor expansion, 
B(x 3) = B—xB,—\3—8) B,+ $x Byp 
+ xd(5 — 8?) B,, + $4°(6-&)P B,, 


* Received in original form 1960 September 5. 
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or, to the second order in x and 6 


178 















B («, 3) =B—xB,—85AB, + $x? B,, + x8A B,, + (AB, + $47 B,,)8, (4) 
where we have denoted B(o, A) by B and defined 


0B 0B 
a-(a)., *(a)., 


with a corresponding notation for the second derivatives. 
Then it follows from (4) that 


if 70) B(x, —. ;) a= B-B,.x—)B,.5+4x°B,, 


+AB,, #5 + (AB, + 327B,,)8 (5) 


where 


\ o(A) B(o, A) da 
tr. [ o(A) B, (0,A) dd 
AB,= |" o(A)AB, (0,A) dd 


B..- | ” ofA) Bez (0,A) dd 
0 





a | ” (A) AB, (0, A) dd 
0 





¥B,, = | * o(A)A2B,, (0, A) dd 
0 
The expression “ for F can be written 
F= {In  o(A) B (#5 5) arm i, o(A) B(o, A) dA—In(1+8)$ (7) 
2E 0 +6 0 


where E = In 10 = 2°303. 
Substituting (5) in this expression we obtain, after some calculations, 


eo? im \S 
P= a {— Bex (1+ Borie Fe - (FB) ] 
2E B B B B 
+ [a+ Be Be (BY |w+ [Be - Babs ao sh. (8) 
7 a (BY 


Substituting this value of F into (1) we have 


logy D=o0'2im—K-—M,-AM-C}+1 .... ~ (9) 
where 
K=K,5+ K,é* 
AM=W,x+W,x* + (10) 


C=C,x8 
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(9) 


(10) 
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and 

K,= 73 (1+ 2B) (11) 

W,= = 7 (13) 

- [7 - Sele. (15) 


The coefficients K, and K, are similar to those given by McVittie (1956, p. 157, 
formula 8.616) but now evaluated at the present epoch. 
K, is in fact Davidson’s «* (1959, p. 62, formula 5.7). 


Using the formula 
- _ 2hP—h, 
x=8| 1 ah? 5... | (16) 


it is possible to write AM as a power series in 5. It is then easily seen that W, is in 
fact Davidson’s A* (1959, p. 61, formula 5.4) and can also be derived from 
McVittie’s expression for AM (1956, p. 157, formula 8.617). 

McVittie’s and Davidson’s treatment of the K-correction and the ‘‘ Stebbins- 
Whitford ’’ effect, which really means any kind of secular variation in the intrinsic 
luminosity of the light source, differ mainly through the way in which the expression 
for F has been split up in an attempt to account separately for the two effects. 

(a) McVittie separates the expression for F into the factors: 


| * (A) B(x, A) dr 


|; o(A) B(o,A) dd 





(17) 





which only depends on the time variable 7, from which the Stebbins-Whitford 
correction is derived, and 


I 7 A 
a [. o(A) B (= 3) dy 











0 (18) 
| o(A) B(x, A) dd 
0 
which reduces to 1 if =o and is used to derive the K-correction. 
(6) Davidson separates the expression for F into the factors: 
I ' A 
a(A) BI ty, .) dx 
5]. ) (m5 (19) 





|; o(A) B(ty, A) dd 


which depends only on 6 and has the advantage of yielding the K-correction 
directly with coefficients evaluated at the present epoch fp, 
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{ R o(A) B (« 4) da 
{ ; o(A)B (‘ 4) dn 


which is only a function of the time if 5 = 0, and from which the Stebbins- Whitford 
correction is derived. 

But (18) and (20) are, in fact, functions of both 6 and x, and if the corrections 
are pushed to a higher order of approximation, one gets also a ‘‘ mixed ’’ correction 
of the form 


and 


(20) 





C=C,xd (21) 
where the coefficient C, is given by formula (15). 

It appears from the above analysis that it might be dangerous to split the 
expression for F from the beginning, since it is impossible, except to the first order 
approximation, to separate completely the corrections due to x and 6. 

So far, only the first order approximations have been used but, if it were 
decided to extend the ‘‘ apparent magnitude-redshift ’’ relation to a higher order 
in 5, a ‘‘ mixed” correction (21) should also be included. 


University of Natal, 
Durban, 
South Africa: 
1960 October 18. 
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